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Abstract 



A number of interesting features of the ground states of quantum spin chains 
are anahzed with the help of a functional integral representation of the sys- 
tem's equilibrium states. Methods of general applicability are introduced 
in the context of the SU(25' -|- l)-invariant quantum spin-*? chains with the 
interaction —P^'^\ where P^^"^ is the projection onto the singlet state of a 
pair of nearest neighbor spins. The phenomena discussed here include: the 
absence of Neel order, the possibility of dimerization, conditions for the ex- 
istence of a spectral gap, and a dichotomy analogous to one found by Affleck 
and Lieb, stating that the systems exhibit either slow decay of correlations or 
translation symmetry breaking. Our representation elucidates the relation, 
evidence for which was found earlier, of the —P^^^ sp'm-S systems with the 
Potts and the Fortuin-Kasteleyn random-cluster models in one more dimen- 
sion. The method reveals the geometric aspects of the listed phenomena, and 
gives a precise sense to a picture of the ground state in which the spins are 
grouped into random clusters of zero total spin. E.g., within such structure 
the dichotomy is implied by a topological argument, and the alternatives 
correspond to whether, or not, the clusters are of finite mean length. 
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1. Introduction 



There is a geometric aspect to the structure of the spin-spin correlations found at low 
temperatures in a number of quantum spin systems. Our purpose is to introduce some 
generally applicable tools for the analysis of such phenomena. That is done in the context 
of the SU(2S' + 1) invariant models introduced by Affleck and also studied by Batche- 
lor and Barber, and Klumper [6,10,11,38,39], which include the spin 1/2 Heisenberg 
antiferromagnet as a special case. 

The systems considered here are one dimensional chains of spin S variables, with the 
Hamiltonian 

H=-{2S+1)Y,J^P^% (1-1) 

X 

where Px^y is the orthogonal projection onto the singlet state of two quantum spins of 
magnitude S, and Jx > 0. The models with translation invariant (all Jx = J) or staggered 
coupling constants (possibly with two different values for x even and odd) are of special 
interest and some of our results are specific for these cases. We will refer to the Hamiltoni- 
ans (1.1) as the spin S models with interaction —P^^\ The explicit form of the interaction 
in the basis of eigenstates of is: 

(25 + l)Pi«)= f2 (-1)"-^IA -/?)(«,-«! (1-2) 

a,l3=-S 

The interaction can of course also be expressed as a polynomial in the Heisenberg interac- 
tion Sa; • Sa;-|-i. For S — 1/2 and S — 1 one obtains 

p(0) _ / 3 - Sa; • Sa;+l for S = 1/2 . . 

^-^y-\l{Sx-Sx+ir-l for 5 = 1 ^^-^^ 
The analogous expressions for general S can be found, e.g. , in [39]. 
The phenomena which we shall address are: 

1 ) The nature of the order parameters which characterize the possible occurence of 
symmetry breaking in the ground state. 

i. Neel order. In higher dimensions such models may exhibit Neel order in the ground 
state. For the standard Heisenberg antiferromagnet this has been proved for dimensions 
d>3 and also for cZ = 2 if 5 > 1 [16,35]. The representation introduced here permits to 
rule out that possibility for the translation invariant models with interaction (1.1) in one 
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dimension, on the basis of known results in percolation theory (in two dimensions). (In 

one dimension, the representation relates Neel order to a transient behavior in a system 
of random loops which form the boundaries of the connected clusters of a random cluster 
model) . 

ii. Dimerization. The one dimensional models may, nevertheless, exhibit a two-fold 
translation symmetry breaking, caused by dimerization. While the interaction favors the 
pairing of neighboring spins into singlet states, not all neighboring spins can be paired 
simultaneously. There are, of course, states - corresponding to different dimerizations of 
the lattice, in which half (or, on the lattice Z*^, a fraction l/{2d) ) of the interaction terms 
are minimized. While these are not true ground states, it turns out that in one dimension 
for S large enough {S > 1) this structure is present in the ground state, which decomposes 
into a superposition of two partially dimerized states. Spins on even sites have stronger 
correlations with their neighbors to the right in one of the states, and to the left in the 
other. 

The classical dimerization picture is too naive in two aspects: 1) the model's cor- 
relation functions extend beyond nearest neighbors, and 2) the spins correlate in larger 
clusters than pairs. A virtue of the method employed here is that it permits to describe 
this phenomenon (and the picture of the correlated clusters) in explicit and precise terms. 
In particular, we find the following behavior: in the state where the spins on the even 
sites are more correlated with their neighbors to the right one finds that with probability 
1 some spins on the left half-infinite chain (— cxd, x] form a singlet with some spins on the 
right half-infinite chain -|- 1, +oo), for each x even. In the same state this probability is 
< 1 for x odd. 

2) The "dimerization versus power law decay" dichotomy. 

The dimerization does not always occur. However, we show that there is a di- 
chotomy: the ground state either dimerizes, or exhibits slow decay of correlations (with 

\x{SqS^)\ — -|-oo). The dichotomy has the following geometric content. When the 
(even) clusters of correlated spins are tightly bound, with only a finite number of clusters 
having the origin in their span, then a topological argument implies that the translation 
symmetry has to be broken. The alternative is that the clusters of correlated spins are 
loosely bound, with the origin (as well as any other lattice site) belonging to the span of 
infinitely many correlated pairs. In that case, the above sum of the correlation function 
diverges. In fact, in our representation, that sum measures the number of correlated spin 
pairs with x < and y > 0. 

The dichotomy discussed here is reminiscent of the one found by Affleck and Lieb for 
the Heisenberg antiferromagnetic spin chain with half-integer spins [7]. However, unlike 
the dichotomy of [3], the one discussed here is not restricted to half- integer spin. The 
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string order parameter mentioned above is also a variant of one which has been found 
relevant before in the context of the Heisenberg model. In fact, the method introduced 
here applies to quite a broad class of antiferromagnetic spin models, to be discussed in a 
subsequent paper [5]. 

3) Relation with the Potts models. 

The ground states of the ID spin chains with the Hamiltonian (1.1), as well as the 
Gibbs states exp{—PH) are related to Potts model, which in case of the translation invari- 
ant interaction are always at the self-dual point. The existence of a relation was noted, 
at the level of a similarity of the spectra of the relevant transfer matrices, by Baxter [8] 
for the spin 1/2 model, and by Affleck (who introduced the general spin S model with 
interaction — Batchelor and Barber, and Kliimper [6,11,39] for general spins. The 
representation employed here makes this relation very explicit. In particular, the dimer- 
ization corresponds to the existence of long range order in the corresponding Potts model, 
and the expectation values of any observable of the spin chain can be expressed in terms 
of quantities calculated within the Potts model. The relation presented here extends also 
to models with inhomogeneous couplings (for which the corresponding Potts models are 
no longer at their transition point), and thus extends beyond the exactly soluble cases 
discussed in [8,6,11,39]. 

4) Decay rate. 

Using the geometric representation, and the FKG inequalities which it allows to bring 
to bear on the problem, we derive an effective bound on the decay of correlations of 
general observables in terms of the truncated two-point function t{x, y) of an associated 
two-dimensional {2S + l)^-state Potts model: 

\{AB^)\<CaCb ^(^'2/) (1-4) 

xGsuppA 
y GsuppBz 

For a more complete statement and the notation see Theorem 7.2. Assuming that the 
truncated two-point function of the two-dimensional Potts model in a magnetically ordered 
pure phase always has exponential decay, this result implies the existence of a spectral gap 
in the ground states of the spin S models with interaction — P*^^) whenever dimerization 
occurs. The case of staggered couplings, { Jeven, </odd}, is of interest for the discussion of 
the spin Peierls instability. Our method confirms the results obtained by Cross and Fisher 
for the exponents describing the leading behavior of the energy and the mass gap as a 

function of | Jeven — JoAa\ [15] . 

This paper serves as an introduction to a technique of wider application, which is based 
on a decomposition of the Gibbs states of a large class of quantum spin Hamiltonians as 
superpositions of what we call quasi-states (see Appendix I). In a subsequent paper we 
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show that such a representation exists for any isotropic nearest neighbour interaction under 
the condition that there are no frustration effects. No frustration essentially means that 
the lattice is bipartite and the Hamiltonian has ferromagnetic interactions only within each 
sublattice, and all interactions between spins of different sublattices are antiferromagnetic. 
In case there are only ferromagnetic interactions the bipartite structure is irrelevant. 

For the more general case, the method used here permits to give a natural description 
of the occurrence of fractional spins at the edges of finite chains. Let us just mention some 
results for the spin-1 antiferromagnetic chain with Hamiltonian: 



with the coupling constants a and h satisfying a > and 6 < 0. Again there is a dichotomy 
[5]: under the assumption of sufficiently fast decay of correlations (expected to be violated 
only when a — —b) either 

1) the ground state of the infinite system dimerizes, and thus breaks the translation in- 
variance of the Hamiltonian, 

or, 2) finite pieces of the chain behave as if near each edge there was an excess spin S = ^. 

The latter case corresponds to the Haldane phase [30] and the spin- 1/2 's at the edges 
have been observed in electron spin resonance experiments on NENP [29,27,28]. 




(1.5) 



X 
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2. Quasi-state decompositions for equilibrium states of quantum spin systems 



In this section we derive the path integral representation for the ground state, and 
the equihbrium states, which is employed in the derivation of the results described in the 
introduction. The discussion applies to a more general class of systems than those covered 
by equation (1.1). We also introduce here the notion of a 'quasi - state decomposition' of 
a quantum state. Some of its basic properties are presented in Appendix I. 

Before turning to the derivation, let us state the net result for the spin S models with 
interaction —P^^\ 

2.1 The functional integral representation for the spin S models with interac- 
tion -P(o) 



Absorbing a convenient constant in its definition, the Hamiltonian is now given by: 



and we are considering finite chains of spins of magnitude S. (The validity of Proposition 
2.1 below is actually not restricted to the one-dimensional case, though it requires the 
model to be frustration free, i.e., to have a bipartite structure). 

We denote by a = {<Tx} a configuration of joint values for the commuting family of 

observables {Sx }• These configurations form a natural parametrization for an orthonor- 
mal basis of the Hilbcrt space of the system's state vectors. We obtain the following 
representation for the matrix elements, in this basis, of the operator e~^^ . 

For each specified pair of configurations, e~^^{a',a) is given by an integral over 
various possible histories of a time dependent configuration a{t) = {ax{t)}, with the time 
t ranging over the interval [0, (3] , and a{0) = a and cr(/3) = a'. What remains to be specified 
is the description of the configurations contributing to this integral, and the measure with 
which they are integrated over. 

The contributing spin configurations are piecewise constant in time. When a change 
occurs, the spins change simultaneously at a pair of neighboring sites, but the two spins 
are constrained to add to zero both before and after the change. A useful description of the 
process is obtained by associating with each spin configuration a collection of 'horizontal' 
bonds (in the space-time diagram in which time is in the vertical direction) linking the 
pairs of related sites in all the discontinuity events. For a technical reason, we find it 




(2.1) 



X 
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convenient to somewhat modify this relation, and consider the bonds as enabhng, rather 
than forcing the spin flips. (This extension yields a higher degree of independence in the 
measure seen below). The integral over the time dependent spin conflgurations is based 
on an integral over those time indexed bonds, which we denote by the symbol ui. 

We shall use the following symbol for the consistency indicator function 



and Iper[c(-)|a;] which equals one if, in addition, a"(-) is periodic in time (a"a;(/3) = (JxiS^))- 

The collection of all the spin configurations which are consistent with a given bond 
configuration uj is conveniently described with the aid of a loop decomposition of the 
space-time diagram (which forms a finite-volume subset of Z x [0,/3] as illustrated in 
Figure 1). The loop to which a point (x,t) belongs is found by moving from it along 
the vertical line till a bond is reached. Upon reaching a bond, the path traverses it, and 
then continues in the reversed time direction along the vertical line to which it just crossed. 
This procedure is continued until the path either closes (by returning to its starting point), 
or reaches the time t = or t = /3. For the time-periodic constraint, the loops do not 
stop at t = 0, /3 but reemerge at the other end. Following these instructions, space-time is 
decomposed into a collection of lines which may form either closed loops or lines connecting 
pairs of 'boundary sites' (at t = 0, /3) (in the non-periodic case). Drawing each bond in 
duplicate, the corresponding lines may be drawn so that they do not cross. For a specified 
a;, the consistent spin configurations are completely characterized by the condition that 
the staggered spins, {—lYaxit) are constant along each of the loops of u. In particular, 
for the periodic constraint, the number of consistent spin configurations is, {2S -\- l)^pe--('^)^ 
with /per('^) being the number of loops. 

The relevant measure for the time dependent spin configurations can be constructed 
by means of a product measure, obtained by the integration over uj with an auxilliary 
Poisson process distribution, p\Q^p]^{duj), and the discrete summation over the (2S' + l)'p'=''*^^^ 
consistent spin configurations (i.e., those with I[cr(.)|ci;] = 1). The Poisson measure is 
characterized by the condition that the mean bond density is 1, and that they occur 
independently in disjoint regions. The contribution of a given bond configuration to the 
integral is enhanced by the factor (25* + l)'pcr('^)^ and therefore its effective weight in the 
partition sum is given by the probability measure 




1 if (Jx{t) satisfies the above described constraints 

and all its discontinuities occur at bonds in u) 
otherwise 



(2.2) 



^ip{duJ) = (^;3)-V[o,^](da;)(25+ l)'-('^) 



(2.3) 
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Figure 1: A space-time configuration oj for the — p(°) quantum spin chain, at an inverse 
temperature (3. As in a more general case, the spins are correlated within loops drawn 
by following the lines in the space-time. A special feature of this interaction is that the 
loops can be viewed as the boundaries of the connected clusters of two random cluster 
models, dual to each other. The shaded areas are the connected A-clusters, the connected 
5-clusters are left blank. Ai, A2, -Bi, -B2 are the four independent regions surrounding a 
bond that appear in the proof of the Euler relation. The trace of the loops on the t = 
line shows a decomposition of the spins into random clusters of zero spin. 
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The situation is summarized in the following proposition. 



Proposition 2.1. For a finite system with the Hamiltonian (2.1): 

{a' I e-^^ \a) = J P[o,/3](da;) ^ l[a{P) = a',a{0) = a] (2.4) 



a{-y.I[a{-)\uj] = l 

a) the partition function is given by: 



Zf3 = Tre^^.(('^+')-^^''^+i-') = j P[o,p]{du;) {2S + l)'p-('-) (2.5) 



Hi) the equihbrium expectation values of observables which are functions of the operators 



S^^ can be expressed as 



i/3 



= J ii{dio)EM) (2.6) 



where fj,{duj) = ^ p[Q^p-^{duj) {2S + l)'p°"'('^) and the expectation functional E^{f) is ob- 
tained by averaging, with equal weights, over all the spin configurations consistent with 



In the above proposition and in the rest of this paper by I[- • •] we denote the in- 
dicator function of the event described between the brackets. Of course, the objects 
Zi3, p[Q p]{duj), E^{f) and //(da;), depend on the size of the finite system and the mag- 
nitude S of the spins. 

The functionals can be extended to the full algebra of observables (see Section 2.5 
for explicit expressions). Thus equation (2.6) is akin to a reperesentation of the equilib- 
rium state as a superposition of states. That, however, is only partially true. The linear 
functionals do not have the full positivity properties of quantum states. Nevertheless, this 
point of view is very useful, and is well justified in so far as the expectation values of the 
special (but important) subalgebra of observables {f{{S^})} is concerned. We refer to 
such functionals as quasi-states. The notion is ellucidated in Appendix I. 

As given by equation (2.7), in each quasi-state E^j the joint distribution of the spins 
takes a very simple form. In particular: 

E^{SlS^) = (-1)1^-^1 C(5)I[ (x, 0) and (y, 0) are on the same loop ] (2.8) 
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with CiS) = ^ EL-5 = ISiS + 1). 

Hence, the spin-spin correlation is proportional to the probability, with respect to the 
effective probability measure on the space of bond configurations, that two sites are on the 
same loop of co: 

{S^Sb = (-l)l^"^IC(<S)Prob^( {x, 0) and {y, 0) are on the same loop) (2.9) 



The rest of this section is devoted to the derivation of Proposition 2.1, and of similar 
results for other systems (e.g. the Heisenberg ferromagnet). The discussion of the specific 
properties of the model with H = — Px^l+i is resumed in Section 3. 

2.2 Poisson process representation of e~^^ 



We now turn to the derivation of the functional integral representation, which is done 
in the broader context of operators of the form 

Hr = -Y,Jbhb (2.10) 

where F is a (finite) collection of sites, i3 is a collection of subsets of F, and for each b E B, 
hf, is a self-adjoint operator acting in the Hilbert space (S)ieb^^' ^^^^ ^« state space 
at the site i, and are non-negative coupling constants. We refer to the sets 6 e B as 
bonds, although for the moment they are not required to be pairs of sites. 



Thermal equilibrium states of the system, and its ground state (approached in the 
limit P oo), are associated with the operator e~^^ . Following is a general expansion 
of such operators by means of integrals over a Poisson process. The symbols p and uj 
appearing here are defined afresh, but their usage is consistent with the example discussed 
in the previous subsection. 

=-'^-^'7*^^^'"^n;,),> (2-11) 

where u> — {(6j, ti)} C i3 x [0, (3] is a configuration of time indexed bonds, J|* is the time 
ordered product: 

TT, hb = hbk - ■■ h^hbi, such that tb^ < ■ ■ ■ < h^, (2.12) 

11 



and p'^ {du}) is a probability measure, under which u forms a Poisson process over B x [0, /?], 
with the Poisson density Y\h Jbdt. Thus, u forms a random countable collection of time- 
indexed bonds which occur independently in disjoint regions of B x [0, /?]. 

The Poisson integral formula (2.11) offers a non-commutative version of the familiar 
power series expansion of the exponential function. 



2.3 The Quasi-state decomposition 



In a wide class of models there exists an orthonormal basis {\a)} of the Hilbert space 
of the system such that for all uj 

(«inL,. hb\a)>0 . (2.13) 

For such models the Poisson integral formula (2.11) provides a starting point for a quasi- 
state decomposition (Q-S decomposition) of the Gibbs state defined by 

W> = ^^ ■ P.14) 

As we shall see, the condition (2.13) can be met, in suitable bases, for both ferromagnetic 
and antiferromagnetic models, and it does not require the existence of a basis in which all 
hb have only non-negative matrix elements as was the case in various previous approaches 
[14,24,26,33,37,44,45,51,52,53] (for a treatment of a much wider class of interactions 
see [5]). 

The Q-S decomposition resulting form (2.13) takes the form: 

(•) = J i,{dw){-)^ (2.15) 

with 
and 

pfo,,l(da;)TVK(c.) 
J P[o,^](d^)T^^(^) 

Important features of the model are reflected in: i) the structure of the quasi-states 
and ii) the relative weights in the decomposition (2.15) of different classes of uj. To 
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illustrate this, we now look in detail at the spin-1/2 ferromagnetic and antiferromagnetic 
models. 



2.4 The spin-1/2 Heisenberg ferro- and antiferromagnet 



The Heisenberg Hamiltonian is 

H^-^ = ^t)lE<,,y>Jx,y^x -^2, + Const. (2.18) 

= - E -^-^vHZ} + (2.19) 
<x,y> 

where a — (cr^, cr^, cr^) are the usual Pauli matrices and the sum is over a set of pairs of 
sites. The signs are chosen such that we can always assume that the coupling constants 
Jx,y are non-negative. The following choice of hi, permits us to cast both the ferromagnetic 
and the antiferromagnetic Hamiltonians in the form of (2.10): 

^ r Tx^y ferromagnet ^ 
^ \ 2Px^y antiferromagnet 

where T^^y and Px^y are the transposition and singlet-projection operators acting in the 
Hilbert space of the sites x and y (so, Tx^y interchanges the states at the sites x and y). 
Use is made here of the relations 

a,-a, = 2T,,,-l = l-4pW (2.21) 



In discussing the matrix elements 

^b\(7') = Yl {^\hbMti){<^tAhb2Wt2)---{<^tk-i\hbkW) (2.22) 



it is convenient to consider a space time picture in which the RHS is viewed as a sum over 
paths in the spin configuration space, with {cr^} defined at all times. That configuration is 
piecewise constant, and the amplitude for the process is determined by the matrix elements 
of the operators {hi,}. Beyond this point the two cases need be discussed seperately. 
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2.4-a: The ferromagnet (H^ ). 



In the ferromagnetic case, the hb are transpositions, which occur with amphtude 1. 
Thus, the time-ordered product (2.12) consists of a sequence of transpositions which result 
in a permutation which we denote 7r(a;). In particular, it is easily seen that 

TrK(a;) = 2'^^'^^ (2.23) 

where ioj) is the number of cycles in the corresponding permutation. (The factor 2 
reflects the dimension of the single-site Hilbert space.) More generally, quantities of the 
form Ti: f{{a^})e~^^ are described by the following construction - which is similar to 
but not quite the same as the one presented in Section 2.1 for the Hamiltonian considered 
there. 

The paths a{t) which contribute are constrained, by the nature of the trace, to have 
a{t = Q) = a{t — 13), and thus are periodic in time. The configuration u is visualised, 
in a space-time graph, by a collection of "horizontal" bonds connecting "vertical" lines 
indexed by the lattice sites. For each a;, the contributing spin configurations are obtained 
by decomposing the graph into a collection of loops. The ferromagnetic loops are obtained 
by replacing the antiferromagnetic bonds in Figure 1 by ferromagnetic ones (see Figure 
2). The loop to which a point (a;, t) belongs is found by moving "upward" along the 
vertical line at that point till a bond is reached. At a ferromagnetic bond h = {x,y}, the 
path crosses from x to y and continues in the positive time direction. This procedure is 
continued until the time t — (3 is reached at which point the path jumps to t = at the 
same site. i.e. t = and t — P are identified. The process is repeated until one comes back 
at the point (a;, t). 



Figure 2: The graphical representation of the ferro- and antiferromagnetic bonds used in 
the drawing of the loop configuration co. 
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The permutation K{u}) and its cycle decomposition is easily read from the above 
picture: K{uj) takes x into the site where the loop drawn starting from (a;, 0) returns, for 
the first time, to t = 0. 



The discussion in Section 2.3 yields now the following proposition. 



Proposition 2.2. For the spin-^ Heisenberg ferromagnet, the partition function and the 
expectation values of observables generated by {a^y^^r are given by: 



2^ 



and 



J = Tre''^<-.^> -^^'y^^^'y-^^ = J pfo^^](da;) 2^^^''^ , (2.24) 



2-F 



j fiFidoj)E^if) (2.25) 



where fi^{du!) is the probabihty measure 

^^^idu;) = (Zj)-Vfo,/3](da;)2^"(-) (2.26) 

and the expectation functional E^{f) is obtained by averaging, with equal weights, over 
all the spins configurations which take common values (r]^ = ±1) on the cycles {7} e u, 
corresponding to K{u}): 

E^{f) = 2-''^-^ J2 mv)) (2.27) 



In fact, E^{f) are quasi-states (adapted to the algebra generated by {al}xer accord- 
ing to the definition in Appendix I). In section 2.5, below, we shall discuss the expectation 
values of other observables. However, let us note here that two interesting choices for 
/ are: /i(cr^) = cr^cTy and /2(o""^) = exp{phJ2xer '^x)- Here, in contrast to our overall 
convention, cr^ denotes a Pauli matrix with eigenvalues ±1. For these functions E^{f) is 
given by 

E^{fl) = l[{x, 0) and (y, 0) are on the same loop] (2.28) 

E^if2) = n cosh(|7|/?/.) (2.29) 
■y€ui 

where I7I denotes the number of times 7 intersects the t — axis (so, I7I = the length 
of the cycle tt^). In fact very similar expressions to the ones above can be derived also 
for itinerant electron models. The analogue of formula (2.29) for the Hubbard model was 
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used by Aizenman and Lieb in [3] to derive a generalization of Nagaoka's Theorem to finite 
temperatures. 

2.4-b: The antiferromagnet (H^^). 

A significant difference between the ferromagnetic interaction and the antiferromag- 
netic one, which is seen aheady at the classical level, is the possibility of "frustration" . Our 
analysis is restricted to the frustration-free case, which is characterized by the existence 
of a bipartite structure: the lattice F decomposes into two sublattices. Fa and Fb, with 
the couplings between two sites restricted to be antiferromagnetic if the sites belong to 
distinct sublattices, and, in more elaborate models in which both kinds of interactions are 
present (see [5]), ferromagnetic within each sublattice. For convenience we also define 

/ 2") 1^-2/1 — /"'"''" ^ belong to the same sublattice ,^ 

^ ' \—l hi X and y belong to distinct sublattices ^ " ' 

As we shall see now, under the assumption of a bipartite structure, the positivity condition 
(2.13) is satisfied even though not all the matrix elements of the operators Kiuj) are 
positive. 

In the computation of quantities of the form Tr /({(T^})e~'^^'*^ we need the matrix 
elements of hb for the antiferromagnet: 

h, = 2P^^^= Yl {-ir-^\P,-P){a,-a\ (2.31) 

a,/3€{|,-|} 

Note the similarity between (2.31) and (1.2). 

At this point we have a choice: the spin 1/2 AF system can be discussed in a form 
close to that of the spin 1/2 ferromagnet or in a form which is suitable for the more general 
spin-S models (with the Hamiltonian (2.1)). In order to lay the grounds for a sequel to 
this paper, where we encounter spin 1/2 systems with mixed F and AF interactions, we 
shall present the first option before treating the general spin case. 

As in the ferromagnetic case, for a given uj the allowed time dependent spin configu- 
rations {cr^} are piecewise constant (in time) and can change only at pairs of sites where 
a bond occurs. When a change occurs it is again a transposition. However, there are the 
following differences from the ferromagnetic case: 

i) there is a restriction that where a bond occurs the two spins add up to zero (both 
before and after the event) 
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ii) at a given bond a transposition may or may not occur 



iii) the amplitude for a given "path" is (-1) raised to the number of transpositions, i.e. it 
equals the parity of the resulting permutation. 

Lemma 2.3. In a bipartite system, for each pair of conGgurations {cr, a'}, all the per- 
mutations which take a into a' and which can be written as products of transpositions 
exchanging sites on different sublattices have a common parity, denoted here sign{a, a'), 
with 

sign{a,a) = +1. (2.32) 



Proof : Each transposition changes the number of positive spins on the A-sublattice 
by ±1. Therefore, noting that the spin configurations take values ±i, the parity of any 
admissible permutation is 

(_l)E.er^(^--<)_ (^2.33) 



We now derive the following formula for the matrix elements of the operators K^uj): 



o-ti.---,o-tfc_i 



= sign(cT, ct')2'o ^""^IafW, a' \ uj] (2.34) 

where Iq^{uj) and IafIct a';uj] are defined graphically, in a decomposition of the space- 
time graph into paths, which in many ways is similar to the one seen in the ferromagnetic 
case. The paths now come in two kinds: open paths with end points of the form (x, 0) 
and (x,/3), and closed paths which we call internal loops (see Figure 1). The open paths 
are constructed by starting from any point (x, t) with x & T and t = or /5, by moving in 
the vertical direction until a bond is met. Upon the traversal of a bond, the orientation of 
the motion in time is reversed. The path stops upon reaching t = or t = p. The internal 
loops are obtained in the same way, starting from any point {x,t) G F x [0,(3] that does 
not already belong to an open path and a loop is completed when the path comes back 
to its starting point. Iq^{uj) denotes the number of internal loops. Iaf[C) cr')"^] is defined 
in terms of the open paths which wc interprete as imposing a pairing condition on the 
spin configurations a = cr(0) and a' = a{P). The condition is that for any pair of points 
X = {x, t) and y = {y, t') with t and t' either or (3, that are the end points of a path in 
a;, one has 

a{t),a{t')y = (-l)l-fl = 1 - 26t,t' (2.35) 
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Iaf[o;o-' I o^] = 1 if (2.35) is satisfied for all open paths in lj and otherwise. 

Formula (2.34) is now a direct consequence of the graphical representation of uj and 
the definitions given above. 

When computing Tr K{uj) one identifies t = and t = (3. Then, all paths are closed, 
i.e. they are loops, and we write l-^^ {uj) or iper(<^) to denote the total number of loops in 

The above discussion leads now directly to the representation given by Proposition 
2.1 for 5=1/2. 

2.5 Structure of the quasi-states 



From the expressions in Proposition 2.1 and 2.2 it is obvious that the quasi-states 
{■)^ depend only on the structural properties of u) revealed in the random loop picture of 
the configuration. In fact, the only relevant property of u is how its random loops link 
together sets of sites at t = 0. In both the ferro- and the antiferromagnet a quasi-state 
is uniquely determined by the permutation ti^ of the sites in F, which takes the site x G F 
to the site ti{x) which is where the loop at x, starting off in the positive time direction, 
intersects t — Q for the next time. There is a one-to-one correspondence between the cycles 
in TTj^ and the loops in uj that intersect the t = hyperplane. Let 70 denote such a generic 
cycle: 

/Xi X2 ■■■ Xr-l Xr\ .2.36) 
\X2 Xz ■ ■ ■ Xr Xi J 

A first observation to make is then that in the functional i?^, with ^ — F or AF, there 
are no correlations between the spins on two sets of sites which support distinct cycles in 
the permutation: 

£^f(nAo)=n^*(Ao) (2.37) 

70 70 

where A^^ is an arbitrary operator acting on the sites {xi, . . . ,Xr}. Moreover, E^{A^q) 
depends on cj only through 70 and we therefore might as well denote it by E^^{Ayg). 

As the E^ are linear, they are completely determined by their values on operators 
of the form A^^^ = cr*^^ ■■■cr*''^, where ij G {0,1,2,3}, and by wc denote the identity 
operator. It is a straightforward excercise to compute E^^^{al\ ■ ■ ■ (j't^^) starting from the 
space-time picture: we consider the spin to be a piecewise constant function along 70 
taking values cr^ = ±1. Where there is a factor the expectation value picks up a factor 
±1 according to what the spin is at that point, reverses the spin and cr^ = ia^a^. For 
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the ferromagnet the result is: 

<« • • • O = Hdk" • • • ^"li) + ■ • -^^1 - h)) (2-38) 

For the antiferromagnet the only difference is that the spins now have a staggered inter- 
pretation: 

al = {-l)\-\a\ al = cj\ al = a' (2.39) 

This is equivalent to reversing the spin whenever the loop 70 changes its orientation with 
respect to the positive time axis, i.e. whenever it traverses an antiferromagnetic bond. The 
relation (2.39) does of course not correspond to a unitary transformation of the full algebra 
of observables of the system and is possible only because we are not dealing with true 
quantum states but only with quasi-states. A compact expression for the antiferromagnetic 
functionals is: 

Ei^{< ■ • •<) = ■ --^t) = (-l)*^"^l'^='^^'^^<« • • •<) (2.40) 

We conclude this section by summarizing the results of Propositions 2.1 and 2.2 in 
the following way: we found a representation of the ferro- and antiferromagnetixc ground 
states as a convex combination of quasi-states which for the ferromagnet are a partition 
of the lattice into subsets on which the spins are locked together in a parallel state. For 
the antiferromagnet the spins are rigidly correlated in a staggered manner. This picture 
can be considered as a generalization of the VBS-states where neigboring spins are paired 
into the singlet state (for a different generalization see [17]). In particular, depending on 
the properties of //(cL;), the states can also have long range order, characterized by the 
fact that the clusters percolate throughout the system. Finally wc want to remark that 
the quasi-states each have a much larger symmetry group (which depends on uj) than the 
ground state itself. Because there are no correlations between the sets of sites belonging 
to different loops, the spins on the distinct cycles in tt^^j can be rotated independently. 

2.6 The SU(2S' -I- l)-invariant spin-s models with interaction — P^^) 

We now turn to the Q-S decomposition of a generalization of the spin-| Heisenberg 
ferro- and antiferromagnet to arbitrary values of the spin. In particular this generalization 
includes the one-dimensional antiferromagnets that are the main subject of this paper. 
Starting from the the Heisenberg Hamiltonians as they were written in Section 2.4, we 
just replace T^^y and 2Pi°y by the corresponding operators for a spin S system: in the 
"ferromagetic" case the interaction = T^^^y interchanges the states at the sites x and y 

and for the antiferromagnet hf^^ = {2S + l)Px'^y, where P^^^ is given in (1.2). It is obvious 
that both and hf^ are SU(2)-invariant. Due to the invariance of this interaction 
under parity preserving relabelings of the 25-1-1 states, this SU(2)-symmetry is actually 
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embedded in a larger (for S' > 1) SU(25' + 1) -symmetry. The model with interaction 
is a degenerate ferromagnet and its ground states are given by all permutation symmetric 
states of the system. (The integrability of the one-dimensional models with interactions 
and ~h[ was demonstrated in the well-known work of Sutherland [50]). The interaction 
h^^ is much more interesting and is the actual subject of the rest of the paper. In the 
case of the one-dimensional lattice, we recover the SU(25' + l)-invariant model that was 
first studied by Affleck [6] and also by Klumper [39] and Batchelor and Barber [11]. 

As before, we can only deal with the non-frustrated case, and let Fa and Tb denote the 
two sublattices defining the bipartite structure of the system. The SU(25' -|- l)-symmetry 
of the Hamiltonian with interaction (1.2) is then represented by the fundamental represen- 
tation on one sublattice and the antifundamental representation on the other sublattice. 

The matrix elements of e~^^ and the partition function for this model can be ex- 
pressed using the same correspondence between configurations u and sets of loops as for 
the spin-1/2 antiferromagnet. But now each loop has to be decorated with a label a taking 
the 25' -I- 1 values —S, —S -|-1,...,<S' — 1,5'. The spin configuration as a function of a; e F 
and t e [0, /?], is now the following: 

^\ _ i OL if the label of the loop is a and x E Ta ^-|^\ 
^ '' ' [ —a if the label of the loop is a and x eTb \ • ) 



The proof of Proposition 2.1 now follows by the same arguments of Section 2.4, taking 
into account the correspondence between loop labelings and space-time spin configurations 
given in (2.41). 

In the one-dimensional case the loops can be interpreted as the boundaries of the 
elements in a partition of the two-dimensional space-time. As we will see in the next 
section, the weights with which these partitions occur arc given by the Gibbs weight for 
the associated configurations of a Potts model at the selfdual point. This will enable us to 
analyse the possible long-range order in the ground state of these models. 
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3. Equivalence with the two-dimensional (25' + l)^-state Potts models 



It will be natural to consider the spin chain with a priori different coupling strengths for 
the even and odd bonds. More generally, we are concerned with the Hamiltonian for a 
spin-S' chain given by: 



where /i^^ = (25* + 1)P(°) is defined in (1.2), > and L_, L+ e Z, L_ < < L+. 

We shall now show that associated to the geometric structure of Section 2 is a Potts 
model, or rather a pair of dual Potts models (the A- and the S-model) . In the translation 
invariant case { J^ = J for all x) one arrives at the Potts model at its self-dual point, where 
it is exactly solvable [8] . In this situation the equivalence was conjectured by Affleck and 
established on the level of the spectrum of the transfer matrices by Batchelor and Barber 
[11] and Kliimper [38,39]. 

The Potts models are defined over a 1+1 dimensional lattice, Z x IR, in which one of 
the directions (corresponding to the "time" of Section 2) is continuous. 

To introduce the lattices on which the A- and 5-Potts model variables reside, we start 
from the space-time of the quantum spin chain embedded in and partition into 
vertical strips of width 1 which we label alternatingly A and B, with the strip < x < 1 
getting the label A, as in Figure 1. The Potts variables of the A-modcl reside on the 
vertical lines bisecting the A-strips and the variables of the S-model are situtated on the 
lines bisecting the S-strips. 

For a volume [0, T] x [L_, L_|_], the Potts configurations of the ^-model are functions 
^(x,t), X = even + ^ G [L_,L_|_], which are piecewise constant in time and take values 
in {1, . . . , g}. In the time direction we always take the periodic boundary conditions. 
Other than that, there are two natural boundary conditions for a Potts model, which are 
exchanged under the standard duality map: the free and the wired, the latter corresponding 
to adding an extra strip to the left and to the right of the volume where the spins are 
rcqiured to assume a common value (we assume the convention that the partition function 
includes also the sum over this common value). As will emerge from our discussion, the 
relevant boundary conditions here depend on the label of the strip along the boundary. If 
it is A then the ^-model gets the free b.c. and the 5-model the wired one, and otherwise 
it is the other way around. 




(3.1) 



x=L- 
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The partition function of the ^-model (with the relevant boundary terms) is given by 



Potts ^ /■^j-(d^)^.^^p^ E / dtJf (%,,,),e(.+2,t) - 1)} (3.2) 

even+1 



where = {J^} and {J^} are sets of positive constants (these are the ferromagnetic 
couphng constants in the vertical and horizontal direction respectively and in which we have 
absorbed the inverse temperature of the Potts model), is a product of independent 
Poisson point processes on the lines {x = 2?i + x [0, T] with intensity J^, and 'Yl^ 
denotes the sum over all configurations ^ for which the discontinuities happen only at 
points (x, t) in the configuration uj of the Poisson process. The sum in (3.2) should be 
interpreted as incorporating the boundary condition convention explained in the previous 
paragraph. 



The Poisson measure incorporates the interaction in the vertical direction, and it can 
be arrived at by way of a continuum limit of ordinary discrete Potts models with vertical 
couplings analogous to those seen in (3.2) for the horizontal interaction. 

For the S-model we have an analogous expression with the condition {x = even +|} 
replaced by {x = odd +|}. 

Potts models are conveniently studied via an embedding in the random cluster model 
formulated by Fortuin and Kasteleyn [19]. As we shall see, it is at that level that the 
correspondence with the quantum spin chain is most explicit. The theorem below and 
Theorem 7.2 present some key results which are derived by that route. 



Theorem 3.1. Let the parameters of the finite quantum spin chain with interaction — p(°) 
on the interval [L_, and the ones of the A- and B-Potts models be related as follows: 



Jf = {2S + 1) J 1 , 4^ = (23 + ir'J^_i 

q={2S + lf,/3 = T (3.3) 



Then 
i) 

-ilnZ(^) = -llnZ---^ Y: J. + 2S J: J. 

L_<x<L_^_-l L_<x<L_^-l 

- ^ (H(-l)''- - (-1)''^} - 1) + 1) (3.4) 
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In particular the ground state energy per site of the quantum chain equals the free 
energy per unit volume of the associated Potts model up to a trivial constant. 

a) The Potts model has periodic boundary conditions in the vertical direction and free 
boundary conditions in the horizontal direction. 

Hi) The distribution //(dw) of the random loop representation for the quantum chain 
as described in Section 2, is identical to the distribution of the boundaries of the 
connected clusters of the Potts model in the FK-representation (see below). 



It should be noted that the couphng constants of both the A and the B Potts model 
depend on all the coupling constants { J^} of the quantum spin chain. The relevant coupling 
constants as defined in Theorem 3.1 are and with x—even +i for A-model, and 
a;=odd +^ for the 5- model. 

In Section 7 we present some explicit relations between the correlation functions of 
the quantum chain and those of the Potts models, which, in particular, imply a non- 
perturbative result on the spectral gap in the quantum spin chains with staggered couplings 
(see Theorem 7.1). 



We now turn to the proof of Theorem 3.1. 



First we derive the Fortuin-Kasteleyn representation [18,19] of the two-dimensional 
Potts models described above, by showing that the partition function and the probability 
measure of the Potts model are equal to the ones of a random cluster model. The random 
cluster model is obtained by considering a configuration uj = u)a x cu^ generated by the 
Poisson process (da;yi)p^ (do's), as a partition of the union of the vertical strips of 
type A in the following way: a point [x^t] G uja is considered as cutting the vertical strip 
with coordinate x at the height t and a point (x, t) E ujb connects the two pieces of A-strip 
adjacent to the S-strip at coordinate x at height t (see Figure 1). The partition function 
of the random cluster (RC) model is then given by: 



= j fji"^ iduA)pi" (duB) g^-(--x--) (3.5) 



where Ca{uja Xi^s) denotes the number of connected clusters (consisting of pieces of strips 
of type A) in the configuration cja x ojb- The equality of Z^°^^^ and can be derived 
straightforwardly by expanding the exponential in (3.2): 

ZT'^ = I pf (da;) 5]-- exp{ ^ I dt jf (5^(,,,),^(,+2,t) - 1)} 
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= J Pa {duJA)pi" {dujB) J2''^ n %2/.t).^(y+2,t) 

€ (y,t)eujB 



,e>),^(x+2,t(r>) 



In the sequel we shall write u for oja x ujb and p'^(da;) or p*^"^ '-^ ^(da;) instead of 
p^''(da;A)p^"(da;B). 

A by now standard and very convenient tool for the study of the randon cluster model 
is provided by the FKG inequalities [20]. Following is the FKG-structure which will be 
used to derive some of the main results of this paper in Sections 4-7. 

We define a partial order on the configurations u as follows: a;' ^ a; if the set of 
A-bonds in u is contained in the set of A-bonds in u' and the set of S-bonds in uj 
contains the set of S-bonds in lv' . (By A-bonds we mean bonds occurring in an A- 
strip and analogously for 5-bonds). As before, let /per('^) denote the number of loops in 
uj, considered with periodic boundary conditions in time. We consider the loops as the 
boundaries of a collection of connected subsets (connected clusters) of the plane. Each 
such connected set consists of the vertical strips n < x < n + 1 connected by horizontal 
bridges. All strips in a given cluster are either of A or of S type. Let Cyi(a;) denote the 
number of connected clusters of A type and Cb{oj) the number of clusters of B type. With 
these definitions one then has the following obvious but crucial properties: Na{oj) and 
Ca{}^) are decreasing functions of u) and Nb{uj) and Cb{(^) are increasing, where Nc{(jj) 
denotes the number of bonds in uj which occur in the strips of type C, C = A, B. 

In the following proposition and throughout the rest of the paper, p'^(da;) denotes 
the Poisson measure on the configurations co for the quantum spin chain on a finite in- 
terval, [L_,L+], containing L+ — L_ + 1 sites, and at inverse temperature (3 (see Section 
2.1). J stands for the collection of coupling constants {Jx}l_<xL+ which determine the 
intensities of the independent Poisson measures for each bond. We will sometimes need to 
distinguish between the coupling constants for even and odd x; we then use the notation: 

J — {Jevenj <^odd)' 

One then has the following relations: 



Proposition 3.2. With boundary conditions described above, the following relations hold: 



24 



for any u > 0, 

= C2p^^~^ (duo) u^^-^Ci^) 

where 

C2 = J{(-^)"--(-^)'+>-^exp{/3((l-«-i) ^ J, + (l-z,) ^ J,)} 

K even a; odd 

C3 = «-i{(-^)'--(-^)"+>-^exp{/3((l-z.) ^ J, + (1-^-1) ^ J,)} (3.8) 

a; even a; odd 



Proof : We start by reformulating the random loop picture as follows. As before we 
consider the loops as being embedded in the plane which we have divided into vertical 
strips labeled alternatingly A and B as in Figure 1. We can then associate with each 
configuration cu, a set of A- and S-clusters as follows. When uj contains an activated bond 
in an A-strip this bond is considered as forming a connection between the two -B-strips to 
the immediate left and right of it and as cutting through the ^-strip to which it belongs, 
and analogously for the bonds in a 5-strip. Thus for each configuration u) we have obtained 
a collection of clusters of A and B type and such that two different clusters with a piece 
of common border arc of different type. C^(a;) and Cb{uj) are the number of A-clusters 
and -B-clusters in lo, respectively. 

We consider a finite interval of the chain of the form [L_,L+], with L_ < < L+. 
and such that the number of sites is even. 

For concreteness, let us suppose first that the number of sites (L_|_ — L_ + 1) is even. 
Then, the boundary strips are either both of the A or both of the B type; suppose they 
are of the A type. Recall that we then have free boundary conditions for the ^-clusters 
and wired boundary conditions for the S-clusters (wired independently at the left and at 
the right boundary) in the space (horizontal) direction. In the time (vertical) direction the 
boundary conditions are periodic. Call everything connected to the leftmost S-cluster the 
"outside" . Then each loop in ut unambiguously encloses a domain and has the complement 
of that domain as its outside, i.e. every loop in uj is the outer boundary of exactly one 
domain and except for the "outside" , each domain has a loop as its outer boundary. 
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Therefore the following relation holds 

IperH = CA{i0) + CB{i0) - 1 (3.9) 

Obviously we would have arrrived at the same relation (3.7) if the boundary strips are of 
the B type or if the number of sites is odd. Denote by Na{(jij) and Nb{(jo) the number of 
activated bonds that occur in to in the A-strips and -B-strips, repectively. Then, with the 
boundary conditions described above, the following Euler relation holds (see Lemma 3.3 
below): 

Cniu) - CAiio) + Na{uj) - Nb{u) = i{(-l)^- - (-1)^+} (3.10) 
Using (3.9) and (3.10) we have: 

= «-i{(-^)'--(-i)'+>-V''(da;) ^2CBM-iVMc.)+iv^(c.) ^g^^^l) 

The factors u^^^^^ and u^^ii^) (.qt^ easily be absorbed in the measure p{duj) as a mere 
modification of the intensity and the normalization of the Poisson process using the rela- 
tion: 

p\duj) A^^'") = e'^^^'--^'^ p^\duj) (3.12) 
for any A > 0. That completes the proof. 



Proof of Theorem 3.1: 

In Proposition 2.1 we obtained the following expression for the partition function of the 
spin S model with interaction — P^°): 

Z^^) = y"p'^(da;)(25+l)'p-('^) (3.13) 

We apply the second equality of Proposition 3.2 with u = 2S -\- 1. Comparing the result 
with (3.6) one then sees that up to a trivial constant Z^*^^ is equal to the partition function 
of the yl-Potts model in the Fortuin-Kasteleyn representation with q— (25 + 1)^. Taking 
the logarithm one obtains i) of Theorem 3.1. In the same way one obtains ii). 



We still have to prove the particular form of the Euler relation that was used in the 
proof of Proposition 3.2. We will use the following lemma also to determine the self-dual 
point of the "continuous time" Potts models. 
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Lemma 3.3 (The Euler relation). With the prescription of above one has: 



E{u;) = Cb{u;) - Ca{u;) + Na{u;) - NbH = |{(-1) 



L_ 



(-1)^+} 



(3.14) 



Proof : It is obvious that (3.14) holds for the configuration u that contains no bonds at 
aU. We now prove that (3.14) is vahd in general by showing that for any configuration u' 
obtained from any other configuration uj by removing a bond, one has E{uj') = E{uj). 

Suppose for concreteness that the bond to be removed is in an A-strip and call it h. 
We divide the neighbourhood of h into four regions and label them Ai,A2, Bi and B2 as 
shown in Figure 1. There are two possibilities: 

1) in the cases where the domains A\ and A2 are connected with the bond h present, we 
must have that after the bond has been removed Bi and B2 are no longer connected. 
So in this case removing b leaves the number of A-clusters unchanged, but increases 
the number of S-clusters by 1. Hence indeed E{u}') = E{u}). 

2) if with h present Ai and A2 are not connected, then there at least one of these two A- 
clusters must be surrounded by a S-cluster that contains both Bi and B2. Therefore 
Bi and B2 will still be connected when h is removed. So, in this case removing h 
leaves the number of S-clusters unchanged and decreases the number of ^-clusters by 
1. Again E{uj') = E{uj). 



The Potts models are exactly solvable at their self-dual point in the thermodynamic 
limit. For quantum spin chains this corresponds to the translation invariant case where 
Jx — J for all X, and taking the limit /3 — > 00. This relation follows from the next 
proposition. 

Proposition 3.4. The relation 



determines the self-dual point of the continuous time Potts model with partition function 




(3.15) 



(3.2). 
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Proof : The statement is a direct consequence of the Euler relation. In Figure 1 a certain 
configuration w of cuts and bonds is depicted. The clusters consist of the pieces of vertical 
strips connected by the horizontal bonds in between. Recall that in the time direction the 
boundary conditions are periodic. In the horizontal direction the clusters at the boundary 
strips are subject to free boundary conditions. We also added an extra strip without bonds 
to the left and to the right of the interval. This is the situation where Lemma 3.3 holds 
and the A and the B Potts model are then exactly each others dual. 

Using Lemma 3.3 and (3.12), the expression (3.6) for the partition function of the 
A Potts model can then be transformed into the partition function of the B Potts model 
with new coupling constants: 

= J p'if{du;A)p%''"{du;B) g^-(--x--) (3.16) 

As (do^yi) generates the horizonal and (da;^) the vertical bonds of the 5-model, 

the condition for self-duality is = qJ^. 



Note that because of the relations between the parameters of the quantum spin chain 
and the Potts models, self-duality of the Potts models is equivalent to translation invariance 
of the spin chain. 

We end this section by deriving an expression for the spin-spin correlation function 
of the spin chain under consideration, in terms of the associated random cluster model 
(3.6). This random cluster model is a simultaneous realization of the A Potts model and 
its dual, the B Potts model. As remarked before, the boundary conditions are periodic in 
the vertical direction for both realizations, but at the left and right edges one has wired 
boundary conditions for one and free boundary conditions for the other realization. 

In terms of the random loop model the (imaginary time) spin-spin correlation function 
is given by (see (2.9)): 

{S^{x,t)S^{y,s)) = (-l)^"^C(5)Prob((a;,t) is connected by a loop to (y,s)) (3.17) 

In the random cluster model the loops are viewed as the boundaries between two domains 
of opposite type. It is therefore clear that: 

{S^ix,t)S'iy,s)) = (-l)--^C(5)//(/A(x,y)/B(x,y)) (3.18) 
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where, for C = A,B, 1^7 (x,y) is the indicator function of the event that and yc 
belong to the same cluster. Here = (xct), with xc defined by the requirements that 
\x — Xc\ = |, and that xc belongs to strip of type C. Note that for any two space-time 
points X and y, y) is an increasing and /b(x, y) a decreasing function. 

We will often use the notation x ~ y to indicate the event that x and y are connected. 
If X and y have integer space coordinates (i.e. they are on the vertical lines) connected 
means connected by a loop or a line in a;. If the space coordinates of x and y are not 
integer the notation means that they belong to the same connected cluster. The event is 
empty if the two points belong to strips of different type or if the space coordinate of one 
of them is an integer and of the other one is not. 

The individual probabilities Prob(xA and y^i belong to the same cluster), can be com- 
puted in the g-state Potts model (see the proof of Proposition 3.2, where we applied 
the Euler relation). So, 

(S'^(x)<S'^(y))(-l)^"^ < C(5')Probpotts(xA and y^ belong to the same cluster) (3.19) 

and if \x — y\ = 1 and t = s, we have 

{S^{x)S^{x + 1)) = -C(S')Probpotts(^ - i and X -|- | belong to the same cluster) (3.20) 



4. Finite systems and the thermodynamic limit 



As a first application of the equivalence with the Potts model obtained in the previous 
section we now prove some preliminary properties of the spin S model with interaction 
_p(o) Except when mentioned otherwise we only assume that the coupling constants are 
strictly positive. 

Theorem 4.1 (Finite systems). For each finite interval [L_, L^], the limit 

0[._,i.|=^ltoOf._..,| (4.1) 

exists, both for the quantum spin chain and for the associated random cluster model. 
Furthermore, if the number of sites (i.e. — L_ + 1) is even, the limiting quantum state 
has Stot — and in the random cluster model there are no infinite lines (almost surely). 



That the finite volume ground state of a finite chain of even length is unique and is 
a spin singlet confirms the antiferromagnetic nature of the models. The above theorem 
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can be viewed as the analogue of the weU known result of Lieb and Mattis [42] for the 
Heisenberg antiferromagnet for the class of Hamiltonians under consideration here. 



Theorem 4.2 (The thermodynamic limit). 

1 ) The finite volume ground states for the quantum spin chains and the probabihty 
measures for the random cluster models converge to a well-defined thermodynamic 
limit provided the parity of the boundary sites is preserved; i.e. for any local observable 
Q the following limits exist: 

(g)even=^ lim hm (g)g. j (4.2) 

1/ _ — >■ — oo , 1/ _ even — L-*-'— 
£/_!_ — t' + oo ,!/_!_ odd 

(g)odd = ^ lim hm (g)g. (4.3) 

Z/ _|_ — »■ + oo , i/ _|_ even 

and with similar limits deGning //even (da;) and //odd (da;). 

2) The relation between the quantum states and the probability measure of the random 
cluster model valid in finite volume (see Sections 2 and 3) persist in the infinite volume 
limit. In particular: 

{Sx ■ Sy)gven(odd) 

(-1)^ yC{S)Prohi^^^^^^^^^^ {{x, 0) is connected by a loop to {y, 0)) 

(4.4) 

3) For translation invariant couplings the states (■)even and (■)odd are translates of each 
other, and each is invariant under translation of period two, as well as under global 
spin rotations. 



Remark 4.3. For translation invariant and staggered couplings the states (■)even and 
(■)odd are also ergodic and weakly mixing [21] under even translations. This fact follows 
from the clustering relation and the clustering properties of the Potts model which are 
implied by FKG arguments (see Section 7 (Theorem 7.2)). 



The rest of the section is devoted to the proofs of Theorems 4.1 and 4.2. 
Proof of Theorem 4.1 

Since the arguments are fairly standard we shall be satisfied with an outline of the proof. 

The existence of the limiting quantum state for the finite chains is trivial. For the 
random cluster model we can use the one-dimensionality of the system and the Perron- 
Frobenius theorem. 
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If the number of sites in the spin chain (L+ — L_ + 1) is even, it follows that if there 
is one infinite line, then there must be at least two infinite lines with opposite parity, i.e. 
such that at any given time the distance between the two lines is odd. In the random 
cluster picture this means that there are two infinite clusters of the same color (A or B) 
which is impossible by the theorem of Burton and Keane [12]. Alternatively, one can show 
by a variational argument that maintaining two infinite lines of opposite parity in a finite 
system costs a non-vanishing amount of free energy per unit of imaginary time. 

By construction the state of the spin chain is necessarily rotation invariant and because 
there are no infinite lines one easily shows that {f{S^oi)) — /(O) for any function /. Hence 
the state must have 5'tot = 



Proof of Theorem 4.2. 

From the FKG structure defined in Section 3 it follows that the measures 
form a monotone sequence in FKG sense (see Appendix II) if L_ and L_|_ are chosen such 
that the boundary strips are always of the same type {A or B). The standard argument 
then guarantees that the thermodynamic limit converges (see e.g. [1]). 

That also for the quantum system the thermodynamic limit converges and satisfies 
the relation (4.4) is not immediately evident because local observables for the quantum 
chain are typically related to probabilities of a non- monotone and non-local event, e.g. 

{S^Sy)even 

= (-1)1^-2^1(7(5') lim Probuji l i 0) is connected by a loop to {y, 0)) (4.5) 

L — * — oo J L even >■ — ' -\-' 

— »- + oo,i/_|_odd 

The event "(a;, 0) is connected by a loop to (y, 0)" is non-local and it is non-monotone 
because 

Probjoops (■'^ ^^'^ y are connected by a loop) 

(X^ and y^^ belong to the same cluster and\ 
^ J ,r u 1 4- 4-u 1 J- ) i^-^) 

and Yb belong to the same cluster J 

which is not the probability of a monotone event since the connectivity of the A clusters 
is decreasing and the connectivity of the B clusters is increasing for the order structure at 
hand. Therefore we need a separate argument to show that 



lim^ Prob^j^ {{x, 0) is connected by a loop to {y, 0)) 

Prob^^^^^ {{x, 0) is connected by a loop to (y, 0)) (4.7) 



L_-^-oc,,L_eveii f'lL_,L 
— ^ + oo,I/_j_odd 
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If (4.7) holds (and also its obvious extensions to probabilities of more general connectivities 
of finite sets of sites), the result follows by Proposition 2.1. 



For convenience we introduce the following shorthand notation: 

L=[L_,L+] 
Am = [-M, M] X [-M, M] 
L_ — > — oo, L_even 
— > +00, L+odd 
Al = Prob^jj^ [x is connected by a loop to y) 

Al,M = Pl'ob^jj^ [x is connected to y in the box Am) 

Aqq = Prob^^^g^ (x is connected by a loop to y) (4-8) 



It is obvious that for each finite L 



Al^ sup Al,m (4.9) 

M 



and also that 



(X and V are connected to Ai r but not^i 
U K 1 • -A A 4.10 

to each other by a loop mside Am ' ' 

For M fixed we take the limsup^^oo of (4.10). For local events the limsup is actually a 
convergent limit. Therefore: 

< limsupyli — Prob^^^g^(a; and y are connected inside Am) 
L— »oo 

< Prob^g^g^ {x and y are connected to A^ but not inside Am) (4-11) 
The RHS is monotone decreasing in M. Taking the limM^o©: 

< limsup — < Vmh^^^^^ix A and yA belong to two distinct OO clusters) (4-12) 
L— »oo 

Since the RHS vanishes, the value of limsup^^g^ is independent of the sequence of L's 
and so the limit exists and equals A^. 
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5. Absence of Neel order 



An infinite volume ground state of a quantum spin chain is said to be Neel ordered if 

liminf(-l)^(So-S^) >0 (5.1) 

r 

The aim of this section is to prove that in the spin S model with interaction — p(^) Neel 
order does not occur. More specificly we prove the following theorem. 

Theorem 5.1 (Absence of Neel order). For translation invariant or staggered coupling 
constants the inGnite volume ground states (•)even and (•)odd satisfy: 

lim (Sn • Sr) = (5.2) 

r^oo^ ^ ^' (odd) ^ ' 



The absence of Neel order in the ground state of isotropic antiferromagnetic quantum 
spin chains is believed to hold quite generally, but a rigorous proof of this general fact is 
lacking. More was done for the the Heisenberg antiferromagnetic chains: 1) for S' = 1/2 
the exact Bethe Ansatz solution shows no Neel order, and 2) for general spin, an argument 
for the absence of Neel order was given in [48] on the basis of a new correlation inequality 
(which is further discussed in [49]). 

A key role for the argument is played by a result of Gandolfi, Kean, and Russo [22] , 
whose adaptation to our system takes the following form. 

Proposition 5.2. For the wired state of the A Potts model 

Prob(within a box of size r, (0, 0) is surrounded by an A-connected path) > 1 (5-3) 

r— »oo 

which holds regardless of q. 



It is easy to see that (5.3) is equivalent to the statement that there is no S-percolation 
in the state obtained with the boundary conditions which favor ^-connections. The reason 
for that fact is that if the A-clustcrs percolate then, then any finite region is encircled by 
A-connected closed paths which pcrvcnt i?-pcrcolation. If there is no A-percolation in 
this A-preferred state then neither do S-clusters percolate. The proof of [22] rests on the 
planar nature of the connectivity graph (i.e. the nearest neighbor nature of the interaction), 
its refiection symmetries, and the FKG property. 
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Proof of Theorem 5.1 



Theorem 4.2 permits us to express the expectation values of any local observable of 
the infinite spin chain in terms of the random cluster measures /Xeven(da;) and //odd(d'^)- 
The correspondence is identical to what Proposition 2.1 provides for finite systems. In 
particular, by (2.9): 

( — l)^(So • Sj-)odd = ProbA-wired((0) 0) and (r, 0) are on the same loop) (5-4) 



As the loops are the boundaries of the connected clusters in the equivalent random 
cluster model, when two sites x and y are on the same loop then both x^i ~ yA (in the 
^-sense) and 'X-b ^ Yb (in the 5-sense). Hence the probability is bounded by 

< 1 — Prob(within a box of size r, (0, 0) is surrounded by an ^-connected path) '>■ (5.5) 

r— >oo 

where the last step is by Proposition 5.2. 



Let us remark that assuming the validity of the exact results for the Potts model, 
absence of Neel order is also implied by the bound on the spin-spin correlation function in 
terms of the truncated two-point function of the Potts model, which is derived in Section 
7. 
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6. Dimerization versus Power Law Decay: a dichotomy 



Despite the result of the last section, the models considered here may exhibit symmetry 
breaking. However, the symmetry is that of translation, and the phenomenon is caused by 
dimerization of the spin chain. The main result of this section is the following dichotomy: 
for the models considered here, the ground state either dimerizes (and exhibits spontaneous 
breaking of the translation invariance), or the spin-spin correlations have a slow (non- 
exponential) decay, satisfying 



Both possibilities occur (see the discussion at the end of Section 7) . 

A version of this dichotomy was first proved by Affleck and Lieb [7], extending a 
result of Lieb, Schulz and Mattis [43]. They prove that the uniqueness of the ground 
state implies the existence of low energy excitations (of the order of the inverse size of 
the system). The argument in [3] applies only to half-integer SU(2)-spin chains, and to 
SU(?i)-spin chains with a self-conjugate representation acting at a site. Thus, the domains 
of applicability of our analysis and that of [7] have some partial overlap, but none includes 
the other. 

We also show that, when translation invariance is spontaneously broken, the two 
periodic states arc distinguished by the nearest neighbour spin-spin correlation function. 
In this sense (and more if one looks into the representation) the states are dimerized. The 
phenomenon can also be detected by the long distance behavior of quantity: 



considered here only for a; < y e Z with x — y odd. 

Observables very similar to Ox,y have been used earlier in studies of ground states of 
quantum spin systems (see [7,25,36,47]) and in the computation of the magnetization of 
the critical Potts model [9] . In the second part of this section we introduce the "total spin 
on half of the infinite chain" as an operator in the GNS Hilbert space of the ground state. 
This quantity appears to us as of more fundamental significance than the string observables 
Ox,y, and the latter can be expressed in terms of it. Moreover this new operator serves 
as a dimerization order parameter which reveals more clearly the detailed nature of the 
dimer order. 




(6.1) 



X 



„■ 27r V'a o3 



(6.2) 
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6.1 The dichotomy 



Following is the main result of this section. 

Theorem 6.1. For the ground states of the translation invariant spin S model with 
interaction —P^'^^ ((3.1) with = J), the following dichotomy holds: 

- either the translation symmetry is spontaneously broken in the infinite volume ground 
states 

- or the spin-spin correlation function decays slowly (non-exponential) with 

J2WlSl)\ = +<x^ (6.3) 

X 

In the Erst case, the symmetry breaking is manifested in the non-invariance of the 
pair correlation: 

(So • Si)even 7^ (Si • S2)even = (Sq • Si)odd (6-4) 

and also in the string order parameter: 

lim (Oo,2JV-l)even = Jim (C>l,2Jv)odd > (6.5) 
iV— >oo iV— >oo 

lim (Co,2Ar-l)odd = Jim (Cl,2iv)even = (6.6) 

A/— >oo N^oo 



Proof : We again rely on Theorem 4.2 to express the infinite volume state in terms of 
the random cluster measures ;Ueven(dc<;) and /Uodd(da;). Since the even and the odd states 
are translates of each other, it suffices to consider one of them, for concreteness say the 
odd state. 

The physical origin of the dichotomy was outlined in the introduction. The random 
cluster representation permits us now to express these ideas in a precise way. 

First we note that the result of [22], on the absence of simultaneous A and B percola- 
tion, implies that the configurations u contain no infinite lines, i.e. the lines (in space-time) 
along which the spins are correlated occur only in the form of finite, non intersecting, closed 
loops. The most important implication is that at any given time, the spins are locked into 
rigidly correlated even clusters, with ^ 5^ = within each cluster. 

The alternative may now be posed as between the following two possibilities: the num- 
ber of loops surrounding each site is either almost surely finite, or it is almost surely infinite. 
Equivalently, there either is percolation or no percolation (and if there is percolation then 
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either A or B percolates, never both, and hence the translation symmetry breaking). The 
zero-one nature of these probabilites is due to the ergodicity of the measures under even 
translations (Remark 4.3). 



In the absence of percolation, by the identity (2.9): 



J2 \{SlSl)\ = C{S) J2 Prob((a;,0)~(y,0)) 




= C{S) Eodd(# connected pairs {{x, 0), {y, 0)}, a; < 0, y > 1) = oo 



(6.7) 



where lEodd is the expectation with respect to the probability measure //odd- 

That the alternative necessitates symmetry breaking can be seen in different ways: 
i) via the existence of either A or B percolation, ii) by the distinction between the two 
sublattices in the values of a topological index, or iii) via the staggered values taken by 
the string order parameter. More explicitly: 

i) If (6.1) fails then the origin is surrounded by a finite number of finite loops. The last 
of those necessarily touches the infinite cluster. For the random cluster state corresponding 
to the spin state (■)odd7 percolation is possible only for the vl-cluster. (The argument is 
presented in Section 5). It then easily follows that 



A relevant question now is whether the probabilities of these events, which are expressed 
purely in terms of the configurations uj, can be expressed in terms of the expectation values 
of some spin observables. As we show below, the answer is -yes! 

ii) An alternative way to express the translation symmetry breaking, though still at 
the level of the random cluster measure, is by considering the topological index: 



The index is defined only if the number of loops is finite, and it is easy to see that, when 
defined, its values alternate (with the overall phase dependent on u>). 

iii) We now show that the string order parameter (6.2), which is an observable of the 
quantum spin chain, directly detects cluster connectivity (as opposed to loop connectivity 
which determines the spin-spin correlation). By an application of (2.6-7): 



Probodd( X belongs to an infinte 




(6.8) 



w{x,u;) = (-1) 



1(7 encloses x) 



(6.9) 




(6.10) 
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with 

EfjjiJDx^y) = I [each loop of 10 intersects [x, y\ an even number of times] . (6.11) 

We use here the observation that the net flux of any simple loop through an interval is 
either zero or ±1. If a; has a loop with a non zero flux through the interval [x, y] then the 
conditional average of Ox y, over the consistent spin conflgurations, vanishes. In the other 
case, EL. ^ 0. 

Hence, the string order parameter can be given a neat geometrical meaning: 

{Ox,y) = Prob(all loops intersect [x, y] an even number of times) 

= Prob(any loop that encloses (x — ^7 0) also encloses (y + ^,0) and vice versa) 

= Prob((a; — ^, 0) and {y + |, 0) belong to the same cluster) (6.12) 

It immediately follows that {Ox,y) = if y — a; — 1 is odd. 



For the odd boundary conditions the percolating cluster can only be of the ^-type 
and it is unique. It follows that 

lim {Ox a;+2Af+l)odd = 1™ Probuo^^((a; — -i, 0) and (a; + 2N + |, 0) belong to the OO cluster) 

_ I for X odd ^3) 
[ for X even 

Equation (6.13) is an explicit proof that the alternative to (6.7) is symmetry breaking. 



Going beyond the last statement, we can also show that the translation symmetry 
breaking is necessarily manifested in the nearest neighbour correlations. By (4.4): 

{SlSl+iU^ = -CiS)PToh,^,,ix 

{J A*odd(dtu)/A(a;, a; +1) if a; is odd 
/ A*odd(da;)/B(a;, a; + 1) if a; is even 

^ f / iiodd{du))lA{x,x + 1) if X is odd , . 

\ / A*even(da;)/A(a;, a; + 1) if a; is even 

where we also used the duality relation. 



When there is percolation, A-clusters percolate in the state //odd but not in the state 
/Ueven- The fact that this difference is then detected also at the level of the nearest - 
neighbor connections is an implication of the general criterion of the Rising Tide Lemma, 
which is derived here in Appendix II. (In the terminology explained there, Ic{x,y) (with 
C = A,B) are strictly monotone functions.) 
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6.2 The dimerization order parameter 



More can be said about the dimerized state in algebraic terms, by refering to the 
Hilbert space associated with the ground states (•)even and (•)odd via the GNS construction. 

In physical terms, when it is correct to view the spins as organized into neutral clusters 
it is natural to talk about the total excess spin to the right of x, i.e., the total spin in half 
of the chain. This observation explains the following claim, which can be derived within 
our representation of the states (•)even and (■)odd- 

Claim 6.2. Under the condition: 

J2\^{SiSl)\ < oo (6.15) 

X 

(i.e., the opposite of (6.1)) the Mowing hmits exist 

Si = hm J2 e"''"""' = lim Slie) , (6.16) 

y>x 

in the sense of strong resolvent convergence of operators in the GNS Hilbert space associ- 
ated with either (•)evenOr (•)odd- 



In terms of the random cluster representation: 

S^ = Sy I[the loop of y intersects (— OO, x] ] (6-17) 

y>x 

We omit here the proof, except for the comment that what is proven explicitely is the 
convergence of the quantities: 

e^^^' ) S)eve„(odd) (6.18) 

for all strictly local spin observables. Equation (6.14) is a natural condition for both 
the existence of the limiting operator S^ and for the proof of the convergence. The first 
statement is the simpler task, e.g., it is easy to see that: 

( \Sl\' )even(odd) < E ■^/('^O'^^') I (6-19) 

y>o 
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The operator is the third component of a vector (under spin rotations), with 
and defined analogously. Using the the invariance of the states: 

( e2S+i ^ )even(odd) = ( -P|5^|2=o /even(odd) (6.20) 

where P^g is the orthogonal projection onto the subspace with l^'^P = 0. 

The operators thus constructed permit to express the string quantity as follows: 

Oa;,y = e^^^--^~^y^ for y> xe Z (6.21) 
Under the condition (6.15), the formula (6.17) for S^, implies the clustering behavior: 

(C'a;,j/)even(odd) ~ {P\s^__^\2=o) {P\Sy\2=o) ~^ (6.22) 

for \x — y\ ^ oo. Hence, by comparison with (6.13) or directly from (6.17), 

( ^|s^|2=o )odd = m , ( ^^15^2=0 )odd = (6.23) 

Thus, the rotation of all the spins to the right of an odd site, by an angle 2k'n /{2S + 
1), A; = 1, . . . , 2S, produces a state orthogonal to the ground state (■)odd- The orthogonality 
expresses the fact that such a rotation will necessarily break an existing bond. That is 
not the case for rotations of the spins to the right of an even site since, within the odd 
state, there is a positive probability that none of the clusters within which the spins are 
correlated are broken by this division. 

Remark 6.3. In the above discussion we used percolation ideas to relate the failure of 
(6.1) to the positivity of m. For half integer spins, that can be replaced by an alternative 
argument based on the operators S^, which shows how remarkable is the fact of their 
existence. For half integer spins 

^-(^isjez) = l[w{x,u;) = +l] = 0,1 (6.24) 

where the important observation is the 0-1 property. Therefore, due to the ergodicity of 
the state //odd (do;) under the even translations, for each x: 

{ )odd = 0, 1 depending on the parity of x (6.25) 

Since, 

Sl = Sl^, + Sl+, , (6.26) 



40 



the projections obey 

^|5JGZ = 1 -^is^+jez (6-27) 
That directly imphes the lack of translation invariance, in the following explicit form: 

(e^''*^-)odd = (i^|s,|ez)odd - (^15^1 ^2+1) odd 

= (-1)^ (6.28) 

where the overall phase was determined by parity considerations. 



This line of reasoning is reminiscent of the structural proof, by Aizenman and Martin 
[4], of symmetry breaking in one dimensional Coulomb systems. It may be noted that a 
string quantity related to the exponent seen in (6.28) was used in the argument of [7], 
which was also restricted to half integer spins. 

The above argument expresses a different mechanism for dimerization than the one 
used in the proof of Theorem 6.1. The restriction to half-integer spin is compensated by 
the fact that (6.28) can be extended to a different class of (frustration free) Hamiltonians 
for which the percolation picture is not valid, where the result does hinge on the parity of 
2S. 
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7. Decay of correlations in the spin S model with interaction — p(°) 

Some elementary properties of the infinite volume limit of the states (■)even and (■)odd were 
discussed in Theorem 4.2. We now present sufficient conditions for the exponential decay 
of correlations in these states. The main result is: 

Theorem 7.1. If either the A or the B Potts model associated with the ground state of 
the Hamiltonian 

H=-Y1 J^Pxl+i^ Jx > for all x, (7.1) 

X 

has an exponentially decaying truncated pair correlation function (of the Potts variables), 
then the limiting states (■)even 3.nd (■)odd exhibit exponential decay of correlations of local 
observables, and a non-vanishing spectral gap above the ground state energy. 



For a more explicit statement we introduce the following terminology. Let Q be a 
local observable of the form 

k 

Q = e-*^ J]5^; e*^ (7.2) 

i=l 

(any local observable can be written as a finite sum of such products). Then suppQ (the 
support of Q) is the set of space-time points {{xi,t)}. We also define suppQ to be the 
minimal interval (consisting of points at time t with integer space coordinate) containing 
suppQ. 

The truncated correlation of two local observables Q and Q , is defined as {Q ;Q ) = 
{QQ ) ~ {Q){Q )■ For the A Potts model the truncated two-point function, in a state 
which is symmetric under permutations of the q values of the spin, has the random cluster 
interpretation 

TA(x,y) = (Cx^CyA) - (CxA)(CyA) 

= c(Q')Prob(x^ and Ya belong to the same connected cluster) C^-^) 

where 

c(g) = ^(g'-i) (7.4) 

with a similar relation holding for the B Potts model. 
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Theorem 7.2. Let (■) denote the expectation in the ground state of a finite chain con- 
taining an even number of sites, or in one of the hmiting states (■)even or (■)odd of Theorem 
4.2. Then, for any pair of local observables of the quantum system of the form (7.2), 

\{Q;Q')\<CqCq: ^ mm^Tc{y,z) (7.5) 

yGsuppQ 
zGsuppQ' 

where Cq and Cqi are invariant under space-time translations of the observables. The 
relations between the coupling constants of the quantum spin chain and the Potts models 
are given in Theorem 3.1. 



The minimum over the two Potts models in (7.5) is important since it expected that in 
any situation one of the correlation functions vanishes, as (x— y) — > oo. If that is indeed the 
case, then the implication is that the limiting states (•)even and (•)odd ^ire always clustering 
and hence pure phases. In the case of translation invariant or staggered couplings we have 
the following remark. 

Remark 7.3. The following inequality is obvious: 

(the connected cluster of X(7 reaches beyondX 
a box of size 2||xc — yc\\ centered at X^ 
(7.6) 

Therefore, whenever one can show that there is no simultaneous percolation of the A- and 
B-clusters, one also obtains 

11 ^™„'^c(x,y) = (7.7) 

In the case of translation invariant or staggered coupling constants the argument of [22] 
applies, but we expect that simultaneous A and B percolation is absent under much weaker 
conditions. 



The bound (7.5), and some of the arguments used in its derivation are reminiscent of 
the estimate derived in [2] for the quantum Ising model in transverse field. However, the 
case considered here is less direct. Prom (3.18) it is clear that the spin-spin correlation 
function of the quantum chain is not equal to the truncated two-point function of the 
associated Potts model. Nevertheless the decay rates of the two are equal (Theorem 7.6). 

At the end of this section we discuss some implications for models with alternating 
coupling constants. 
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There are three steps to the proof of Theorem 7.2: 

1) When the truncated correlations of quantum spins are transcribed in terms of the 
Potts model, we obtain two distinct contributions. The first is easily bounded in terms of 
the random cluster model's connectivity function, and the second is a truncated correlation 
function of suitable observables of the random cluster model. 

2) Using a general domination principle for FKG measures, the latter correlation 
function is bound in terms of four point functions of the form 

(/c,(x,y);/c.(u,v)) (7.8) 

where Cj are either A or B, and 7^7 (x,y) are indicator functions for the events that xc 
and yc are in the same connected C-cluster (for the definition of xc see at the end of 
Section 3). 

3) It is shown that for any combination of Ci and C2 the truncated correlation (7.8) 
is bounded in terms of the connectivity function of the ^-model. By symmetry, the same 
is true for the S-model (in any situation, only one of these bounds will have a non-trivial 
consequence). The argument is based on a combination of the FKG inequality and the 
Markov property of the random cluster measure. 

The argument relies on two auxilliary results. The first one, used in step 2, is the 
following domination principle for correlation functions with respect to an FKG measure. 

Lemma 7.4 ([46]). For a pair of functions f and g, let F and G be two monotone functions 
with which f + F and g + G are increasing and f — F and g — G are decreasing. Then the 
truncated correlations with respect to any measure with the FKG property satisfy: 

\{f;g)\<{F;G) (7.9) 



The proof is elementary: it consists of two linear combinations of the four correlation 
inequalities resulting from the monotonicity of / ± F and g ±G. 

The next Lemma is needed for step 3. It provides an upper bound for the truncated 
correlation of two monotone increasing functions of the random cluster model. 

Lemma 7.5. Let Ii and I2 he two monotone increasing events for the random cluster 
model, which are determined by the A-connected clusters of two non-random sets Di and 
D2. Then: 

< (/i ; h) < Prob(Di and D2 are ^-connected) (7-10) 
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Proof : The key observation is the foUowing inequahty which results from the combination 
of the Markov structure and the FKG property of the random cluster measure: 



E[/2|/[2?17^D2]/l]<(/2) (7.11) 

The reason for this inequality is that in the complement of Ca[Di], the ^-connected cluster 

of Di, the system can be considered as having free boundary conditions on the boundary. 
These boundary conditions mask the positive event Ii which is determined within Ca[-Di]. 

The above reasoning was first employed in [13] for the Ising model and similar argu- 
ments have later been used for a variety of other applications. 

Using (7.11) and the trivial bound < 1, one then obtains: 

= E[/i I h] = I /[Di ~ L»2]/2]E[/[Di ~ D2] I I2] 

+ IE[/i I I[Di rf. D2\h\W\Dx rf. D-2\ I h\ 
< E[/[L>i ~ D2] I h] + ih) (7.12) 



Thus, 



(/i)(/2) < (/1/2) < (/i)(/2) + E[/pi 7^ D2]] (7.13) 



where the first inequality is just by FKG. That implies the bound stated in the Lemma 
for the truncated correlation (/i;/2). 



Proof of Theorem 7.2: 

For convenience, we first carry the analysis for observables which are products, of the form 
(7.2), of only and variables at non-coincidental sites, and break the proof into the 
steps described above. 

1) In terms of the random cluster representation of the spin chains: 
{Q;Q') = J ^i{du)E^{QQ') - (^j ii{du)E^{Q)^ (^j ^idw')E^,{Q')^ . (7.14) 



Therefore 



\{Q;Q')\ < 



+ 



fi{du;)E^{QQ') - j f,{duj)EUQ)EM') 
J fx{du;)EUQ)E^{Q') - Q f,{du;)EUQ)^ Q ii{du;')E^>{Q') 



(7.15) 
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The first term on the right side can be interpreted as the average over oj of the truncated 
correlation within the quasi states E^. The second term is a truncated correlation function 
of the random cluster model. 

To estimate the first term we note the factorization property (see (2.37)) 

E^{QQ') = E^{Q)E^{Q') if lo does not contain a loop (7.16) 

connecting suppQ and suppQ' 

Using also the bound \E^{Q)\ < ^I^^pp^I, one has 



li{du;)E^{QQ') - J i^{dio)EUQ)E^{Q') 

< ^IsuppQI + lsuppQ'l J2 Prob(x and y are connected by a loop) (7.17) 



xGsuppQ 
yGsuppQ' 



When X and y are connected by a loop then both ~ y^ and x^ ~ ys, where xc ~ yc 
is our notation for the event that xc and yc belong to the same C cluster. Thus we can 
continue the bound as: 



< ^IsuppQi+isuppg'i y- Prob(xc~yc) 

xGsuppQ 
yGsuppQ' 

<C(?)-^ ^IsuppQI + lsuppQ'l ^ j^^j^ Tc(y,z) (7.18) 



C=A,B 
yGsuppQ 

zGsuppQ' 



2) In order to estimate the second term in the RHS of (7.15) we invoke Lemma 7.4. 
For the functions 

f{uj)=E^{Q), g{u)^E^{q') , (7.19) 



we take 



, x,yesuppQ 



and 



G(a;) = 25l^-PPQ'l ^ /^(x, y) - /^(x, y) . (7.21) 

yx,yGsuppQ' J 

That the conditions of Lemma 7.4 are satisfied follows from the following observations: An 
elementary change of a configuration oj consists of the addition or the removal of a single 
bond. Such a change will either join a pair of loops or cut a loop into two. If E^{Q) is 
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affected, the change is by not more than 2S^^^^^^^. However, the change is zero unless there 
is a pair of sites x, y G suppQ whose loops are either joined or disconnected in the process. 
In this situation the value of the increasing function yesuppQ -^-aI^) x) ~ y) is 

changed by at least 1. 

Applying the Lemma we get: 

J tJi{duj)E^{Q)EM) - M(da;)£;,(Q)^ (^ j y.{duj)E^{Q')^ 
< 45|suppQ|+|suppQ'| ^ ^ |(/c,(x,y);/c,(u,v))| (7.22) 

C\ ,C2 x,yGsuppQ 
u, vGsuppQ' 

which concludess the second step of the proof. 

3) Our goal now is to estimate the quantities (/ci (x, y); /c2(u, v)) appearing in the 
right side of (7.22). By the "A-S" symmetry, it suffices to derive an estimate in terms of 
the ^-model's connectivity function ta(-, •)• 

The events /^(x, y) are of two possible types: 

i) two A-sites, x^ and ya, are in the same connected ^d-cluster, or 

ii) two -B-sites, x^ and y_B, are not separated by a connected A-cluster (i.e. are 
connected by a S-cluster). 

In either case, the events are determinded by the A cluster of the set [xc, yc]^ defined 

by: 

[xa, Ya] a = the two point set {xa, Ya} (7.23) 
[xb, Yb] a = the collection of the A sites in the interval joining x^ and ys (7.24) 



For the application of Lemma 7.5 we define two monotone increasing events as follows: 
if Ci = ^ put Ii = /A(x,y) and Di = [xA,yA]A; if Ci = 5 put = 1 - /A(x,y) and 
Di = [x.B,yB]A- h and D2 are defined in terms of /c2(u, v) and [x(72,yc2]A in the same 
way. Lemma 7.5 then implies: 

|(-fci(x,y);/c2(u,v))| < Prob([xci,uci]A ~ [ucj, vc2]a) 

< J2 Prob(zc~z'c) (7.25) 

where we have also used that | (1 — Ii ; 12) | = | (-^i ; -^2) | 
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Interchanging the roles of A and B and using (7.3) one finally obtains the estimate 
|(/c,(x,y);/c,(u,v))|<c(g)-i^mm^ J] tc(z,z') (7.26) 



Combining the inequalities (7.15), (7.18), (7.22), and (7.26), one obtains the estimate 
(7.5), for the case where Q and Q' are products of and operators at distinct sites. 

For the general case we use the identity S'^ = i[S^,S^] to express the product of 
spin operators as a linear combination of products of only and S^. The products may 
contain repeated factors. With a trivial modification the argument given above applies to 
such products as well. 



Proof of Theorem 7.1: 

Theorem 7.1 is a direct consequence of Theorem 7.2. 



It is interesting to consider the implications of the above analysis for the spin chains 
with alternating coupling constants 

J J J even If 3^ cvcn 

^ \ <^dd if X is odd 

First it should be appreciated that while the spin system's Hamiltonian is only periodic, 
the associated Potts models arc translation invariant. A great deal is known about such 
Potts models (though at varying levels of mathematical rigor), and our methods allow one 
to extract from that some relevant information on the ground states of the quantum spin 
chains. 

For the Potts models it is known that the correlation length is a meaningful notion, 
defined by the limit 

ePotts = - Alnrc(0,a;) (7.28) 

The existence of the limit is an easy consequence of the supermultiplicativity (itself implied 
by an FKG type argument): To(a;, z) > to{x, |/)to(|/, z), where tq — c{q)~^Tc, C — A or B. 

Since Theorem 7.2 provides only upper bounds, let us first strengthen the relationship 
between the long distance behavior of the correlations in the two models. 
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Theorem 7.6. For the models with alternating coupling constants (7.27): if ^Potts < 
then also 

^QSChain = " ^ In | (^o^^,^) | (7.29) 

exists and ^Potts = ^QSChain- 



Proof : Let < m — ^Potts defined by the limit (7.28). The theorem wiU follow from 
the following inequalities: for any e > not too large, there exist constants C+,C_ > 
such that 

^^^_rn{l-e)\x\ > | ^^3^ ^3^ | > ^_^-m{l+e)\x\ ^7 3(3) 

for all X large enough. Since the spin-spin correlations are dominated by the function 
t(-), the upper bound is trivial (in fact, an auxilliary argument shows that C(5')e~'^l^l 
would do). We now provide the argument for the lower bound on the spin-spin correlation 
function. For concreteness let us assume that there is no percolation in the A Potts model. 

By (2.8) the problem amounts to estimating from below 

P(0, x) = Prob(0 and X are on the same loop) 

We first show that for all small e > and finite D > {D will be taken 0(1)) the following 
quantity satisfies an exponential lower bound as in (7.30) 



Pe,D{^,x)=PTOh 



/there are sites fX G [(— 4£a;) a, 0] , f G [x, {x + 4:ex)A],\ 
such that U and V are on the same loop and each of 
the ^-strips at the edges of the two intervals does not 

^contain any bonds at times t G D] j 



(7.31) 



with X odd. As the outer boundary of any A-cluster is a loop for the spin chain we 
obviously have: 



^^£,^(0,2;) >Prob 



/there is an ^-cluster intersecting both [— 4£a;, 0] and [x, X -\-\ 
Aex] but not [—AeXjX + AexY, and each of the ^-strips 
centered at {—Aex)A^ — |, and {x + 4:ex)A does not 

ycontain any bonds at times t G [—D, D] j 



and therefore 

(3 an ^-cluster connecting [— OO, 0] and [x, -|-Oo] and each of^ 
the ^-strips centered at {-^Ex) a-, \-,^-\-, and {x^A.ex)A 
does not contain any bonds at times t G [— -D, Z)] 

— Prob(3 an ^-cluster connecting [— OO, 0] and [a; -|- 4£a;, +Oo]) 

— Prob(3 an ^-cluster connecting [— OO, — 4£a;] and [a;, +Oo]) 



(7.32) 
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i) By FKG and translation invariance 

(3 an A-cluster connecting [ — OO, 0] and [x, +Oo] and each of\ 
the A-strips centered at (— 4ex)/i, |, and {x+4:ex)A I 

does not contain any bonds at times t G [—D, D] ' 

> Prob(3 an ^-cluster connecting [—OO, 0] and [x, +Oo]) 

X Prob(the yl-strip centered at ^ does not contain any bonds at times t G [~D, -D])^ 
and, using (7.28), 

Prob(3 A-cluster connectmg [-00, 0] and [x, +Oo]) > c(q)" Va(0, x) > C' e~^^+^^'^^''^ (7.33) 

ii) The middle term in (7.32) can be replaced by the straightforward estimate 

Prob(3 A-cluster intersecting both [— OO, 0] and [x + 4ex, +Oo]) 

< c"e-(^+2=)^l^l (7.34) 

where C" > depends only on m, and it is assumed that < £ < |. By symmetry, (7.34) 
applies also to the last term in (7.32). Combining (7.32-33) we obtain the desired bound 
for Pe^£)(0,x): for some C > (depending on D): 

PsM^^^) > C(l - Ce-=^l^l)e-^^+^)^l^l 

In order to complete the proof we have to argue that the exponential lower bound on 
Pe,D{0,x) implies a similar lower bound on P{0,x). More precisely, we use the bound for 
-?£,£) (0, x), which we only need for "0" even and x odd, to obtain the desired estimate on 
P{y, z) with y = — 1 or y = and z = x or z = x + 1. This covers all possible combinations 
of parities for the two sites. 

We use Lemma 7.7 given below to show that there is some o: > such that for y = — 1 
or and 2; = a; or a; + 1: 

Prob(?/ and Z are on the same loop | the event described in (7.31)) > e""^"^'^' (7.35) 

from which the desired bound immediately follows. 



The basis for the claim (7.35) are the following three observations: 
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i) The event described in (7.31) implies the existence of a hne (forming part of a 
loop) outside of the rectangular neighbourhoods Bi = [{—4ex)A,0] x [—D,D] and B2 = 
[x, {x + 4ex)^] X [— -D, -D], connecting a point of the form {u, s) with a point {v, t), where 
u e [{—Aex)A-i 0], f G [x, {x + 4:ex)A] and s^t — ±D. 

ii) For every pair of points (u, s) and {v, t), on the boundary of the boxes Bi and B2 
respectively (see i)), there is a pair of local events Si and S2 in the boxes Bi and B2 such 
that under Si (w, s) is connected by a line to (—1,0) or (0,0), ad libitum, and under S2 
{v,t) is connected with {x,0) or {x + 1,0). The events, whose choice depends on u,s,v, 
and t, are depicted in Figure 3. 

iii) The conditional probability of the above mentioned local events Si and S2 occurring 
together - conditioned on any explicit configuration in the complement of Bi and B2 - is 
not smaller than e""^"^!^', for some a > 0. 

More explicitely, the local events are constructed as follows (see Figure 3). Assuming 
the line connecting Bi and B2 (in the complement of bi and B2), reaches Bi at the point 
Uo{u,D). The event ^^i(uo) inside Bi is then described by specifying that in the vertical 
strips between u and one sees alternatingly the following picture: in the first, third, fifth, 
... strip, counted starting from u, there is no bond at times < t < D and at least one 
bond at a time —D < t < 0; in the second, fourth, ... strip there is at least one bond with 
< t < D and no bond with —D < t < 0. The strip to the left of u is required not to 
contain any bonds at times —D<t<D. It is then obvious that in the box Bi there will 
be a line connecting {u,D) and (—1,0) and (0,0). The case m = is treated by a trivial 
modification of the above prescription. The event ^2(vo) in B2 is defined in a similar way. 



Figure 3: A schematic representation of the local event Si used in the proof of Theorem 
7.6. Crossed areas do not contain any bonds and gray areas contain at least one bond. 
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Note that the random loop measure jJL conditioned on an arbitrary configuration out- 
side any finite domain in space-time is of the form stated in Lemma 7.7. The function 
/ is the number of loops inside the finite volume taking into account the connections in 
the configuration outside, and q — 2S + 1. It is then obvious that / always satisfies the 
bounded-gradient condition with a = h = 1 because the addition or removal of a single 
bond can change the number of loops by at most one. 

The indicator function of the event 8i fl 82 is of the form FG with F increasing and 
G decreasing. 

It is useful to introduce the auxiliary events ^uo.vq; which form a partition of the event 
described in (7.31): {uq = (tt, s),Vo = {v-it)} is the "first" pair of points on the top or 
bottom boundaries of Bi and B2 connected by a line outside these boxes ( "first" e.g. in 
lexicographic order). We also define Jmin(max) = niin(max){ Jeven, Jodd} and put D = 1. 

The observations i-iii) and Lemma 7.7 then yield for y — (—1, 0) or (0, 0) and z — (x, 0) 
or (a;+ 1,0): 

Prob(y and Z are on the same loop | the event described in (7.31)) 
> ^ Prob(£:i(uo)n£:2(vo) I ^:uo,vo)Prob(£:uo,vo I 

the event described in (7.31)) 



uoeesi 

V0 6SB2 



2[e\xW 



> jy p'^™/^^'^+^^(da;)I[ there is at least one bond in the strip [0, 1] X [0, 1]]| 
{ / there is no bond in the strip [0, 1] X [0, 1]]| 

I j p^^'^+^^'^'"""(da;)I[ there is no bond in the strip [0, 1] X [-1, 1]] | 

l)Jmax ^(^l _ g-Jmin/(2S + l)^g-(2S + l)J„,ax|^^^''^'^ 



2[e\x\] 



X 



X 



^ g-4(2S+ 



for some a > independent of e and x. We can now conclude that for £ > small enough 
one has 

for all X large enough, which is a lower bound of the form (7.30). 



In the above argument we use the following Lemma, in which we refer to the order 
structure on the space of configurations given by the inclusion relation: a; > a;' iff the set 
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of bonds in to contains the set of bonds in uj'. The proof is a rather standard FKG-type 
argument, and is therefore omitted here. 



Lemma 7.7. Let //(da;) be a probability measure of the form: 



//(do;) 



where pj^{du) is a Poisson measure on configurations u in a Gnite volume A (in our case, 
A is a subset of the of space-time), q > I, and fiuj) is a function of bounded gradient, 
in the sense that there are constants a,b > such that f + aN is non-decreasing and 
f — bN is non-increasing with N{lij) the total number of time-indexed bonds in uj. Then, 
for any two functions, F non-decreasing and G non-increasing, which depend on u in 
two disjoint subsets of A (one determining F and the other determining G) one has the 
following comparison inequalities, with expectation values with respect to modified Poisson 
measures: 

[ pidiv)Fiiv)Giiv) < [pi'idiv)Fiiv) [pl''idiv)Giiv) 

(7.36) 

J p{du;)F{oj)G{u;)> j p^' {du)F{u) j pi\du)G{u) 



Assuming now the vahdity of all the results on Potts models presented in references 
[31,40,41,54] (not all of which have been derived rigorously), we obtain the following 
implications. 

i) If Jeven 7^ JoAdi the ground state is unique, with exponential deacy of correlations and 
a spectral gap. 

A perturbative version of this statement, for small (or large) enough ratio of the two 
couplings, is contained in a theorem of Kennedy and Tasaki [36]. 

ii) The case Jeven = Jodd corresponds to self dual Potts models, in which case there is 
a dichotomy, which we have discussed above in Section 6. Its manifestation in the 
Potts model language is: at the self dual point there either is a first order phase 
transition with coexistence of the ordered and disordered phases, or the transition is 
second order - with a unique state at which the correlations decay by a power law. 
The threshold value of g is g = 4, which corrsponds to = 1/2. The implication of 
the Potts model calculation is: 

ii.a For S = 1/2, there is a unique state, but with algebraic decay of correlations 
and no spectral gap (equation (6.3) is satisfied). For the spin system that was 
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known independently of the Potts model result, by the Bethe Ansatz solution 
and a result of [43] . 

ii.b For S > 1/2 the translation invariance is spontaneaously broken and there are 
two partially dimerized ground states, translates of each other, with exponential 
decay of correlations and a non-vanishing spectral gap (the exact value of the 
spectral gap can be calculated, see [9,38,39]). 

iii) In the case 5' = |, for weakly staggered coupling constants, Jeven = 1 + 5, JqAA = 
1 — 5, < 5 « 1, the correlation length is finite but a divergent function of 5. 
The conjectured values of the corresponding exponents of the 4-state Potts model 
[54], yield the following behaviour for the ground state energy per site e and the 
correlation length ^ 

£-£0~-|(^|3, C~l<^r3 (7.37) 

(up to logarithmic corrections). 

This behavior of the spin | Heisenberg antifcrromagnetic chain with alternating cou- 
pling strengths was first obtained by Cross and Fisher in their study of the spin-Peierls 
transition [15]. It has the implication that when the ellastic deformations of the underly- 
ing lattice are taken into account, the ground state of spin-| Heisenberg antifcrromagnetic 
chain develops the spin-Peierls instability. (For a related rigorous result see e.g. [34]). 
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Appendix I. Quasi-state decomposition for quantum states 



Often, basic properties of the state of a quantum system are elucidated by presenting 
it as a convex combination of states with a particularly simple structure. In this work we 
find it useful to consider a broader class of affine decompositions - into convex combinations 
of what is called below quasi-states. These are linear functionals which are required to 
meet the positivity requirements (which are part of the definition of a quantum state) only 
in their restrictions to certain Abelian subalgebras. 

Quasi-state decompositions made already an implicit appearance in the discussion of 
the itinerant ferromagnetism in reference [3]. The utility of such decompositions there, 
stems from the fact that the ferromagnetic (or antiferromagnetic) properties of a system 
with rotation-invariant spin-spin couplings can be expressed through the correlation func- 
tions of a commuting family of the spin variables (e.g., {cr^} with x the site index and 
3 referring to the third component of the Pauli spin matrices). The structure of the re- 
striction of the state to such families is made particularly transparant by decomposing it 
into a convex combination of states in which the correlations are either "0 or 1", in the 
sense that to spins are either locked in a parallel (or antiparallel) state, or are completely 
independent. This is also the major characteristic of the Fortuin-Kasteleyn representation 
of the Ising model [19,18]. 

The purpose of this section is to formulate this notion, and discuss some of its general 
properties in the context of a simple example. 

La Definition of quasi-states 



Definition I.l. Let A be a sub-algehra of obscrvables of a quantum system. A densely 
defined linear functional p (providing the expectation value for observables) is called a 
quasi-state relative to A, and we will say that it is well adapted to A, if p is: 

i) normalized: p(I) — 1 (I.l) 

a) positive on A, i.e. for each bounded A e A, p{A*A) > (1.2) 



The property of states which is not required of quasi-states is the general positivity. 
I.e., the unrestricted validity of (1.2). 
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For finite dimensional quantum systems, the observables form matrix algebras, and 
the linear functionals take the form: 

p{A)^TtQA (1.3) 

where Q is a corresponding density matrix. When A is the maximal- Abelian algebra of 
matrices which are diagonal in a certain common basis {Icin)}, the conditions on Q under 
which p is a quasi-state relative to A are: i) TrQ = 1, and ii) Q has non- negative diagonal 
elements, 

(ttn I <5 I ttn) > forn = 1, ... (1.4) 

The last requirement is considerably weaker than the condition ii') Q > (as a matrix) 
needed for p to be a state. 

Remark 1.2. At this point it is natural to ask whether (1.4) has any A— independent 
content. The answer is rather negative: any hermitian matrix with TrQ>0 satisGes (1.4) 
in some basis. 



As we shall see, in certain situations a convex decomposition of a state into quasi- 
states (relative to a naturally relevant algebra), provides a great deal of insight into a 
state's structure. Such decompositions take the form: 

P = ^PnQn (1-5) 

n 

with the Qn quasi states relative to a common algebra , and pn weights satisfying: 

Pn>0, 5^Pn = 1 (1-6) 

n 



I.b. An elementary example 



As the simplest demonstration of the notion introduced above, consider the system 
consisting of two spin-1/2 objects, with the spin-spin ferromagnetic interaction: 

H=-ai- a2 (1.7) 

The ground state of the related antiferromagnetic Hamiltonian {—H) is given by the rank- 
one projection 

Po = |0)(0| (1.8) 
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onto the singlet state (|0) = (||,— |) — | — |, |))/-\/2), where the total spin is S' = 0. In 
this state one can safely say that Bi = — ^2, e.g. , 

TrPo(^i +^2)^ = (1.9) 

The ground state of the ferromagnet is slightly less elementary. It has the three-fold 
degeneracy of the space on which S =1. The corresponding projection is Pi = I — Pq = -S' 
(S takes here only the values and 1), and the (normalized) state operator is 

Q+ = \S (I.IO) 

While the spins are as parallel as the uncertainty relations allow, it is not true that Bi —B2 = 
0. In fact, 

Trg+(ffi -ff2)^ = 1 > (I.ll) 

The fact that except for the "zero-point fluctuations" the spins are basically aligned, is 
easily seen in the following quasi-state decomposition of the ferromagnetic ground state: 

Q+ = |(i)I+i(iT) (1.12) 

where T = 25 — 1 is the exchange operator (T|a, h) = |6, a); a, 6 = ±). 

Two basic observation here are: 

- For any non-zero vector n G IR"^ both 1 1 and \T are quasi-state operators relative to 
the Abelian algebra .An generated by {Bi • n, ^2 • n}. The former is actually a state 
operator. 

- On any of the Abelian factors An, the quasi-states corresponding to |l and to \T 
represent two simple alternatives. In the first state, Bi ■ n and B2 ■ n form a pair of 
uncorrelated variables (taking the values ±^ independently of each other), whereas 
in the quasi-state |T, Bi ■ n and B2 ■ n are locked together. 

By the rotational symmetry, it suffices to derive these statements for n = (0,0,1), 
within the standard Pauli spin matrix representation, which makes the first statement 
really obvious. For the second we note that the restriction of Tr(iT)A to A G A(o,o,i) is 
determined by only the diagonal elements of in the basis which diagonalizes {o"^,^!}. 
That diagonal part of \T is \{\\, \){\, \\ + \ — \, ~\){~\^ ~\\)- particular, it follows 
that the quasi-state \T yields the following ferromagnetic analog of (1.9): 

TriT(ai -^2)2 =0 (1.13) 

Thus, insofar as the restriction to An is concerned, the ground state of the quantum 
ferromagnetic Hamiltonian is equivalent to a positive temperature state of a classical ferro- 
magnet, and the decomposition (1.12) is identical to the Fortuin-Kasteleyn random cluster 
decomposition of the latter. 
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I.e. Words of caution 



Some notes of caution are due here. The lack of full positivity implies that quasi-states 
may lack some of the properties which are familiar for the expectation value functionals 
(•) associated with regular states, such as: 

i) A>0^{A)>0 (1.14) 

ii) (A) < \\A\\ (the supremum norm of A) (1-15) 

iii) {A*A) = 0^ (AB) = for any (bounded) B (1.16) 
In fact, the quasi-state operator |T provides us with the following "counter-examples" : 

i) Tr ir[(ai - ^2)^ - 1] = -1, while (^i - ^2)^ - 1 > (1.17) 

ii) Tr |T[2 - (ai - ^2)^] = 2, while ||2 - (^i - 0=2)^ || = 1, (1-18) 
and 

iii) TV lT{af - air = 0, yet Tr lT{af - ai){T\ - i, i|}[= 1] ^ 0. (1.19) 

These elementary assertions follow from (1.13) and the observation that 2— (ffi— 0^2)^ = 
2S-1 = ±1. 

Nevertheless, the restriction of a quasi-state to any of the algebras for which it is well 
adapted is free of the above "pathologies". In particular, the three principles (1.14)- (1. 16) 
are satisfied as long as A and B are restricted to any common 

One may also note that while a quasi-state is not fully characterized by its restriction 
to a single Abelian factor, in the above case the operator T is uniquely determined by its 
rotation invar iance and the values of Tr ^TA for A e 

Appendix II: Strictly monotone observables for Gibbs-fields and the Rising 
Tide Lemma 

The FKG structure of the probability measures /Ueven(da;) and /Uodd(dc<;) is quite es- 
sential for our analysis in Sections 4-7. Here we provide some details concerning this 
structure. We also introduce the notion of weak strict monotonicity and prove the Ris- 
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ing Tide Lemma which was needed in Theorem 6.1 to show that breaking of translation 
invariance {peven 7^ A^odd) imphes staggered behaviour of the nearest neighbor spin-spin 

correlation ((Sa; • Sa;+i)even(odd) (Sa;+1 " Sa;+2)even(odd))- 

Recall from Section 3 the definition of the partial order on the set of configurations 
uj = (a;^,a;B), with liJa{b) the set of ^(5)-bonds labeled by space-time points. The partial 
order is defined by declaring 

UJ if loa 3 ^'a UJb C ijj'b (II-l) 

That leads to the notions monotonicity of functions and domination of measures: i) a 
function / is called monotone increasing if 

j{uj) < f{oj') for all u^u' (II.2) 

ii) a probability measure //(do;) is said to be dominated by another probabiltiy measure 
^'(da;), denoted as //(da;) ^'(da;), if 

J //(da;) /(a;) < J ^'(da;) /(a;) for all monotone increasing / (II-3) 



This structure is useful in a number of ways. 

1) Fortuin, Kasteleyn and Ginibre [20] provide a general criterion, under which a 
probability measure would have the property: 



/ 



//(da;) /(a;)^(a;) > y //(do;) /(a;) J f^idco) g{u;) (II.4) 



for all monotone increasing / and g. The FKG condition is satisfied by the random cluster 
measures (see [1]). 

2) For /i satisfying (11.4), the probability measures conditioned on the configuration 
in the complement of a finite volume, are the monotone in the sense: 

E^(- I <^A=) < 1E^(- I ^Ac) for all uja'^ ^ (H-S) 



3) There exist domination relations among the random cluster measures corresponding 
to Potts models with various values of q {q=l corresponding to independent percolation). 
Applications of these relations can be found in [19,20,18] and [1]. 



4) Of particular interest in this work is the domination relation 

A^even 

(da;) ^ //odd (da;) 
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(11.6) 



which follows by a simple limiting argument from (II. 5). 

The last relation includes the statement (Sq ■ Si)odd < (So ■ Si)even- Our goal now 
is to explain why whenever the two states are different, the inequality holds in the strict 
sense. A key role in the argument is played by Lemma II. 2, whose name draws on the 
observation that a rising tide lifts all the ships in the harbor. 

It is useful to introduce the notion of strict monotonicity . 

Definition II. 1. Let Ai be a set of prohahility measures on a conGguration space O with 
an order structure as above. A monotone increasing function on O is strictly increasing in 
the weak sense with respect to Ai (for short M.-strictly increasing), if 

J IJ,{duj) f{Lo) < J v{du)) f{uj) for all fiji' e M, fJ, ^ I' andji^v. (II-7) 



The difference between the two random cluster measures is due only to the boundary 

conditions (pushed to infinity) . They share, however, a common rule for the finite volume 
conditional distributions. (A phenomenon exhibited also by the family of the Gibbs equi- 
librium states at a first order phase transition.) In the terminology discussed in [23] the 
measures have the same set of specifications. 

Lemma II. 2 (The Rising Tide Lemma). Let A4. be a family of probability measures 
with common specifications 7 = {1E(- | a;A<:)}A- A sufficient condition for a function f to 
be Ai - strictly monotone increasing is that beyond some finite volume Aq the conditional 
expectation is strictly monotone increasing with respect to the boundary conditions: 

1E(/ I UAc) < nf I ^Ac) (11.8) 

for all A D Aq and all pairs of boundary conditions such that 

<^Ac <^Ac and uj\c ^ uj'f^c mod 7 , (II-9) 

where the last inequality means that the induced conditional expectations are not identical. 



Remark II. 3. For a given set of specifications, the relevant notion of boundary conditions 
consists of the equivalence classes of {u>ao} defined by 

ujao=uJao iff 1E(- I cuAc) = 1E(- I a;;^c) (ILIO) . 
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It is easy to see that for the random cluster model the equivalence class of boundary 
configurations is determined by specifying the connectivity relations of the boundary sites 
via the connecting paths in A'^. (Obviously many conhgurations in the complement of A 
would be equivalent.) 



Proof of the Lemma: Let and f e M., n ^ i^- Then, by the general result of HoUey 
([32]), there is a coupling p which is a probability measure on the space of pairs {(a;, a;')} 
with marginals /i and v which is supported on pairs with to z^^' ■ Using first the conditional 
expectation formula and then the coupling, we have: 



Assuming the measures are difi'erent, there is some finite volume A D Aq for which p 
assigns a positive measure to pairs of configurations which are inequivalent as boundary 
conditions in A^. Since the integrand in (11.11) is strictly positive at such points, and is 
non-negative in general, it follows that z^(/) — //(/) > 0. 



Lemma ILS. For the random cluster model, for any x, y e Z x H, the random variable 
/yi(x, y) satisfy the criterion of the Rising Tide Lemma II. 2. (The set Aq can be taken as 
any box with x and y in its interior.) 



Proof : Let A be a finite box containing x and y in its interior, and let ujac ^ i^\c- Then 
each of these configurations induces a partition of the A- sites on the boundary of A into 
clusters connected in the exterior. The relation between the configurations means that 
the boundary partition corresponding to uj\c is a refinement of that corresponding to . 
The difference in the conditional expectations of //i(x, y) is the result of two effects: i) the 
induced measures in A are different, and ii) the finite volume conditions under which x and 
y are connected are different (since u' provides a better connected boundary). Denoting 
by I^^^c and 1^' ^ the corresponding indicator functions (of only u!\), we have: 



E(/^(x,y) I c^;.) - lE(/^(x,y) | cja^) 

= ]E(I,^^ (x, y) I o^^c) - E(I,^, (x, y) | uja^) 

= ]E(I^^, (x, y) I uj'j^c) - E(I^^, (x, y) | uja^) + E(I^^^ (x, y) - I^^, (x, y) | c^^c) 




> El(I(x ~ y w.r.t. a;^c and x 7^ y w.r.t. lo\c) \ a;^c) 



(IL12) 
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When the two configurations are not equivalent then, due to the strict positivity of the 
measure, the last term is strictly positive. 



In this paper the Rising Tide Lemma was applied in Section 6. Wc expect it to be of 
interst also in various other situations where the FKG structure is of relevance. 
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